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Abstract / Zusammenfassung 



Balloon releases are one of the main attractions of many fairs. Helium filled rubber balloons 
are released to carry postcards over preferably long distances. Although such balloons have been 
considered in atmospheric sciences and air safety analysis, there is only scarce literature available 
on the subject. This work intends to close this gap by providing a comprehensive theoretical 
overview and a thorough analysis of real-life data. All relevant physical properties of a rubber 
balloon are carefully modelled and supplemented by weather observations to form a self-contained 
trajectory simulation tool. The analysis of diverse balloon releases provided detailed insight into 
the flight dynamics and potential optimisations. Helium balloons are found to reach routinely 
altitudes above 10 km. Under optimal conditions, they could stay more than 24 hours airborne 
while reaching flight distances close to 3000 km. However, external weather effects reduce the 
typical lifetime to 2-5 hours. 



Ballonwettbewerbe sind eine der Hauptattraktionen von vielen Volksfesten. Heliumgefullte Gum- 
miballons werden mit ciner Postkarte gestartet, um eine moglichst grofie Strecke zuriickzulegen. 
Obwohl solchc Ballons in der Atmospharenforschung und der Luftfahrtsicherheit eine gewisse Rolle 
gespiclt haben, gibt es nur sparlich Literatur zu diesem Themenkreis. Diese Arbeit versucht diese 
Liicke mit einem umfassenden theoretischen Uberblick und einer sorgfaltigen Analyse von realcn 
Mefidatcn zu schliefjen. Alle relevanten physikalischen Eigenschaften eines Gummiballons werden 
sorgfaltig modelliert und mit Wetterbeobachtungen erganzt, um ein abgeschlossenes Programm 
zur Trajektoriensimulation zu erstellen. Die Analyse verschiedener Ballonwettbewerbe eroffnete 
detaillierte Einsichten in die Flugdynamik und mogliche Optimierungen. Es zeigt sich, dafi Heli- 
umballons regelmaBig Hohen um 10 km erreichen. Unter optimalen Umstanden bleiben sie mehr 
als 24 Stunden in der Luft und konnen dabei fast 3000 km weit fliegen. Jedoch beschranken externe 
Wettereinfliisse die typische Lebenszeit auf 2-5 Stunden. 
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Kapitel 1 

Einleitung 



Ballonwettbewerbe sind ein beliebter Programmpunkt auf vielen Feiern und Volksfesten. Ein He- 
liumballon wird mit einer Karte und einer Riicksendeadresse gestartet, urn eine moglichst grofie 
Entfernung zuriickzulegen. Der Erfolg ist dabei nicht gewifi - zu unberechenbar erscheint das Wet- 
tcrgeschehen, und zu gutcr Letzt muB die Karte auch noch gefunden und zuriickgeschickt werden! 
Was aber passiert mit dem Ballon nach dem Start in der Luft? Welchen Einflussen ist er un- 
terworfen, und lafit sich eine gute Platzierung mit einer grofien Flugweite vielleicht doch geziclt 
beeinflussen? 

Es lohnt sich, mit einer weiter gefassten Perspektive zu beginnen. Ballons haben schon friih eine 
entscheidende Rolle bei der Erforschung der Atmosphare gespielt [T3] . Zuerst wurden Heifiluft- und 
Gasballons eingesetzt, um Forschern die direkte Untersuchung der Atmosphare zu ermoglichen. 
Spater riickten vermehrt unbemannte Ballons in den Vordergrund, die eine Vielzahl automatischer 
Messungen bis hin zu astronomischen Beobachtungen durchfiihren konnten. Mit automatischen 
Mefiballons (sogenannten Radiosonden) gewonnene Wetterdaten bilden noch heute das Riickgrat 
der taglichen Wettervorhersage |14j . Im Zuge der Entwicklung immer ausgefeilterer numerischer 
Wcttcrmodcllc und Simulationen zur Verbreitung von Luftschadstoffen werden Ballons auch zur 
„Markierung" von Luftmassen verwendet, um die Richtigkeit von Stromungsmodellen |741 [3] zu 
iiberpriifen. Die Umsetzung kann auf mehrere Weisen erfolgen: Durch Ballons, die automatisch 
einem konstanten Druckniveau folgen 54, 5Tj . durch die Auswertung von regelrechten Ballonwett- 
fahrten [7] bis hin zur Einbindung von Ballonwettbewerben in die Forschung 91 . Eine ausfuhrliche 
Ubersicht iiber den Stand der Trajektorienberechnung ist in |92j enthalten. 

Die Untersuchung der Dynamik von Kartenballon^] spielt meist nur eine untergeordnete Rolle. 
Spezielle Arbeiten wurden zur Umweltvertraglichkeit von Ballons [15| und der maximal erreich- 
baren Hohe durchgefiihrt }52) . wobei aber nicht nach einer umfassende Beschreibung eines frei 
flicgcnden Ballons gesucht wurde. Die Eigenschaften des Ballongummis sowie der absolut geringe 
Auftrieb eines Kartenballons erschweren die einfache Ubertragung der charakteristischen Eigen- 
schaften grofierer Ballons. Eine Ausnahme bildet die Arbeit von Roberts [52], die ein einfaches 
numerisches Modell zur Flugweitenabschatzung vorstellt. Ein Vergleich mit realen Kartenballon- 
fliigen oder die Berechnung von Ballontrajektorien erfolgte in dieser Arbeit aber nicht. 

Die offenen Fragen sind also noch zahlreich: Wie weit kann ein Kartenballon fliegen? Wie lan- 
ge bleibt er in der Luft? Wie hoch fliegt er? Kann die Flugweite optimiert werden, und wenn 
ja, wie? Diesen Fragen soil in der folgenden Arbeit systematisch nachgegangen werden, um ein 
umfassendes Bild von dem Verhalten von Kartenballons in der Atmosphare zu zeichnen. Zuerst 
erfolgt eine detaillierte Untersuchung aller relevanten physikalischen Eigenschaften eines Karten- 
ballons in Kapitel [21 gefolgt von einer Zusammenfassung der wichtigsten Umweltcinfliissen in den 
Kapiteln [3] und 2J Diese Bausteine werden anschliefiend in Kapitel [5] zu einem Computermodell 
zusammengefugt, das dann mit Hilfe von Ballonwettbewerben und Experimenten in Kapitel [S] 
gepriift wird. 

1 Im Folgenden werden die bei Ballonwettbewerben verwendeten Heliumballons kurz als „ Kartenballon" bezeich- 
net. 
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Kapitel 2 

Ballonphysik 



2.1 Einleitung 

Die theoretische Beschreibung eines Kartenballons umfaBt alle fiir die Trajektorienberechnung 
relevanten Eigenschaften: Angefangen von der Geometrie des Ballons, dem Materialverhalten des 
Ballongummis, iiber Auftrieb und Aerodynamik bis hin zu dem scheinbar unmerklichen Verlust des 
Fiillgases wahrend des Fluges. Fiir die meisten dieser Fragestellungen ist es ausreichend, alleine 
den Gummiballon zu betrachten. Die angebundene Karte verandert neben dem Gesamtgewicht 
nur die Aerodynamik des Kartenballons. Erst nach dem moglichen Platzen des Ballons wird das 
Geschehen ausschliefilich durch die frei fallende Karte bestimmt. 

Die folgenden Abschnitte gehen detailliert auf die Eigenschaften eines Gummiballon ein, wobei 
die Aerodynamik unter Beriicksichtigung der angehangten Karte in Abschnitt 12.61 diskutiert wird. 

2.2 Geometrie 

Die Grundlage fiir die Beschreibung der Ballongeometrie ist die Kenntnis des Volumens V und der 
Oberflache A. Ein Ballon kann grob durch einen prolaten Ellipsoid mit langer Halbachse a und 
Durchmesser D beschrieben werden (siehe Abb. l2~Tj) . Wird die Abweichung von diesem Ellipsoid 
in einem Parameter e zusammengefafit, ergibt sich das Volumen zu: 

47T 

V = —sR 2 a, R = D/2 (2.1) 

Bci gleich gefertigten Ballons ist das Verhaltnis a/D in etwa konstant und die obige Formel kann 
vereinfacht werden zu: 

V = fvD 3 (2.2) 
A = f A D 2 (2.3) 

Mit Hilfe von Ballonphotos wurden die freien Parameter fiir zwei unterschiedliche Fabrikate be- 
stimmt: 



Mittelung iiber 


- fv fA 1 - R/a 


48 melloc© 32 cm Ballons 
10 Volksbank Werbeballons 


0,96 ± 0,02 0,63 ± 0,02 3, 71 ± 0, 12 0,20 ±0,03 
0,94 ±0,01 0,58 ±0,01 3, 48 ± 0,05 0,15 ±0,02 



Die beiden Fabrikate unterscheiden sich hauptsachlich in der Abplattung 1 — R/a, die in der 
letzten Spalte mit angegeben ist. Die Fehler geben jeweils die Streuung innerhalb der Ballonstich- 
proben an. Fiir manche Anwendungen ist es niitzlicher, das Verhaltnis von Oberflache zu Volumen 
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Abbildung 2.1: Ballon mit dem Aquivalenzellipsoid mit Durchmesser D und langer Halbachse a. 

anzugeben. Von alien Korpern hat die Kugel das kleinste Verhaltnis: 
A/V 2/3 = 4,84 Kugel 

A/V 2 I Z = 5, 05 ±0,20 melloc© 32 cm Ballon (2.4) 
A/V 2 I Z = 5, 00 ±0,09 Volksbank Ballon 

Einen ersten Hinweis auf die elastischen Eigenschaften eines Ballons gibt das Fiillvolumen, das zum 
Platzen der Gummihiille fiihrt. Um eine erste Ubersicht zu erhalten, wurde eine klcinc Stichprobe 
ausgewahlt und bis zum Platzen gefiillt. Wahrend des Fiillvorgangs wurden die geometrischen 
Abmessungen regelmafiig erfafit. Die daraus abgeleiteten Platzvolumina sind: 



Dmax[cm] 


^max [1] 


Vo[l] 


Bemcrkung 


34 ±1 


23 ±2 


0,1 


Volksbankballon bei 20 °C 


37 ±1 


32 ±3 


0,15 


melloc© 32 cm Ballon bei 20 °C 


34 ±1 


25 ±2 


0,15 


melloc© 32 cm Ballon bei -10°C 



Alle Ballons haben sich kurz vor dem Platzen teilweise erheblich deformiert, so dafi die einfache 
Anwendung der Gleichungen 12 . 1 1 und 12 . 2l nicht zuverlassig ist. Hinzu kommt, dafi die Versuchsbe- 
dingungen eine nur ungenaue Annaherung an die langsame Ausdehnung eines steigenden Ballons 
sind. Hystereseeffekte in der Gummihiille und insbesondere die Wechselwirkung mit den fallenden 
Temperaturen beim Aufstieg wurden nicht erfafit. Aus dem Volumen Vq der ungcdchnten (d. h. 
druckfreien) Ballons kann die maximalen Dehnung des Ballongummis ermittelt werden. Sie liegt 
im Bcrcich von 500 % bis 600 % (bezogen auf den ungedehnten Ballon) und ist mit den Angaben 
in |15] vergleichbar. 



2.3 Innendruck 

Die Gummihiille eines Ballons ubt einen zusatzlichen Druck auf das Fiillgas aus. Die Beschreibung 
der nichtlinearen elastischen Eigenschaften von Gummi (siehe z. B. [53]) erlaubt die Berechnung 
des Drucksprungs Ap als Funktion des Ballonradius R. Das Ergebnis der von [64] durchgefiihrten 
Rechnung ist 

Ap(R) = *Y" Cl (^-^)(l + ^^) (2.5) 



3 V a \R R 7 J V Ci Rl 
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Abbildung 2.2: Balloninnendruck in Abhangigkeit des Ballondurchmessers. 



Vn = 



47T 



(2.6) 



mit dem Schubmodul C\ « 0,17 MPa und dem nichtlinearen Parameter C~\jC\ sa 0, 1. Rq ist der 
Radius des entspannten (d. h. druckfreien) Ballons. Das Volumen Vcummi der Gummihiille kann 
iiber die Dichte des verwendeten Gummis aus dem Ballongewicht abgeleitet werden: 



^Gummi — PGummi^Gummi 

Typische Gewichte fiir die bereits genannten Fabrikate sind: 



(2.7) 



^Gummi 


melloc© 32 cm Ballon 
melloc© 13 cm Ballon 
Volksbank Werbeballon 


5,1±0,2 g 
0,86 g 
3, 2 ±0,3 g 



Die Sticliprobe der 13 cm Ballons war zu klein, urn eine representative Streuung zu ermitteln. Das 
Minimum und Maximum der Druckkurve (siehe Abb. 12. 2[) ist fiir C-\jC\ = 0, 1: 



• max = 1,476 xRq 
r min = 3, 143 x i? 



A/w = 0, 745 x 
Ap min = 0,632 x 



4% 



3 V 

u m mi 

3 



Vn 



(2.8) 
(2.9) 



Mit der bekannten Abhangigkeit des Innendrucks von der Ballongrofie kann nun das Volumen 
eines mit der Gasmenge N gefiillte Ballons bei einem vorgegebenen Druck pg und Temperatur To 
berechnet werden. Unter der Verwendung der Zustandsgleichung fiir ideale Gase erhalt man fiir 
die Bcstimmung des Ballonvolumens V die Gleichung: 



(po + A P (V))V = NRT 



(2.10) 



R ist hier die universelle Gaskonstante. Obwohl der Drucksprung (Gl. I2.5|) selbst keine monoton 
steigende Funktion des Ballonvolumens ist, trifft dies auf die gesamte linke Seite zu: 



— (( Po + Ap(V))V) >0 V V >0 



(2.11) 
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Abbildung 2.3: Vcrcinfachtes Strahlungsgleichgewicht, bestimmt durch die Sonneneinstrahlung J 7 , 
die diffuse Strahlung .Fir + cos(#)<xF von der Erdoberflache und die Warmestrahlung des Ballons. 



Dies garantiert fur alle realisierbaren Kombinationen der freien Parameter jeweils eine eindeutige 
Losung fiir das Ballonvolumen. Da der Druck iiber einen weiten Bereich nur wenig variiert, wird 
er in dieser Arbeit durch einen konstanten Wert angenahert: 



Ap 



3 V 



Ci 



4760 Pa x 



^Gummi 

/3 cm 3 
Fo/100 cm 3 



(2.12) 
(2.13) 



2.4 Thermodynamik 

Die Temperatur eines Ballons bestimmt das Volumen iiber die Zustandsgleichung des Fiillgases 
und die elastischen Eigenschaften der Gummihiille. Wahrend die Temperatur eines Ballons im 
Schatten in guter Naherung durch die Lufttemperatur gegeben ist, mufi in der freien Atmosphare 
und insbesondere in grofien Hohen das Strahlungsgleichgewicht mit der Umgebung berucksichtigt 
werden. 

Der Arbeit von Nelson 68 folgend wird eine einfache Zweibandbeschreibung fAbb. l2~3|) gewahlt. 
Im Bereich optischer Wellenlangen erhalt der Ballon Energie durch die direkte Strahlung der Son- 
ne T und den durch die Albedo a der Erdoberflache reflektierten Anteil. Im Bereich infraroter 
Wellenlangen strahlt der Ballon selbst Warmeenergie ab und empfangt Warmestrahlung J^r von 
der Erdoberflache. Die umgebende Luft der Temperatur To tragt iiber direkte Warmelcitung bei. 
In der folgenden Bilanz wird die adiabatische Abkiihlung bzw. Erwarmung des Ballons beim Auf- 
steigen und Sinken vernachlassigt: 

a 0p t7ri? 2 (l + 2 cos(6»)a)J r + a m 2nR 2 J" IR = 2vrA L uftNui?(T - T ) + 4Tra m R 2 aT i (2.14) 

Die Sonneneinstrahlung erfolgt unter dem Zenitwinkel 9, T ist die zu bestimmende Temperatur 
des Ballons und R sein Radius, a ist die Stefan-Boltzmann-Konstante. Im Folgenden wird ein 
graues Ballonmodell mit Absorptionskocfhzicntcn aiu = 0, 5 und ao P t = 0, 25 als plausible Werte 

1 Die Solarkonstante ist 1367 W/m 2 , Weltorganisation fiir Meteorologie (1982) 
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angenommen. Realistische Werte der Albedo a liegen im Bereich 0, 2 ... 0, 5. ALuft ist die moleku- 
lare Warmeleitfahigkeit der Luft. Die wirkliche Warmeleitung ist um den FaktoiH Nu grofier. Eine 
empirische Formel fur die Nusselt-Zahl einer Kugel ist nach [51] : 

Nu = 2+(0,4Re tot 1 / 2 + 0,06Re tot 2 / 3 )Pr ' 4 f^iy /4 (2.15) 

Der Unterschied der dynamischen Viskositat der Luft im freien Raum rj^ zu dem Wert auf der 
Kugeloberfiache rjx ist fiir die hier betrachtete Anwendung so gering, dafi dieser Korrekturfaktor 
vernachlafiigt werden kann. Die totale Reynolds-Zahl Retot setzt sich aus der Umstromung der 
Kugel und dem konvektiven Anteil, beschrieben durch die Grashof-Zahl Gr, zusammen: 

Re tot 2 = Re 2 + 0,4Gr (2.16) 

Gr _ jr-^g (217) 

In der letzten Gleichung ist g die Erdbeschleunigung, v die kinematische Viskositat der umgeben- 
den Luft und D der Durchmesser der Kugel. Bei einer Anstromung mit 1 m/s erhalt man fur einen 
Ballon etwa Nu ~ 100. Fiir eine einfachere Behandlung wird der Strahlungsflufi aus der Umgebung 
in Gl. 12.141 durch eine effektive Infrarot-Temperatur T e g- ausgedriickt: 

4a m aT c 4 ff = 2A Luft Nu/.R(T - T ) + Aa IK aT 4 (2.18) 
« 2A Luf tNu/i 1 , (r-To)+4ai R (7r o 4 + 16 a i R aT 3 (T-r ) (2.19) 

Umstellen der Warmebilanz ergibt: 

T_T ° " NuA Luft /(2 ff ( ra IR i)+4T 3 (2 ' 20) 

Fiir die betrachteten Temperaturdifferenzen T — To konncn die Strahlungsterme um To entwickelt 
werden: 

T — T ~ ^ cS — — C9 9FI 

10 ~ NuA Luft /(8T 3 aai R i?) + 1 K ^ L) 

Die Temperaturdiffercnz T e g — To reicht von Werten um — 40°C (Nachts, grofie Hohen) bis liber 
+40°C (Bodennahe, Sonnenschein) . Um den Einflufi auf die Ballontemperatur abzuschatzen, wird 
der Faktor in Gl. 12.211 ausgewertet. Die Warmeleitfahigkeit ALuft ist unabhangig vom Druck und 
nur proportional zur Wurzel der Temperatur, so dafi die Warmeleitfahigkeit am Boden auch fiir 
grofiere Hohen ein guter Richtwert ist. Mit R = 15 cm ergibt sich die Abschatzung: 

NuALuft = 2>8x Nu x /2MKy (2 . 22) 

Der durch die Ballongrofie bedingte starke Einflufi der dircktcn Warmeleitung rcduzicrt den Tem- 
peraturunterschied zur umgebenden Luft auf Werte um 10 K. Fiir die Berechnung des Ballonvolu- 
mcns und des Auftriebs ergeben sich nur Korrekturen im Prozcntbereich, so dafi die Ballontempe- 
ratur T gleich der Temperatur der umgebenden Luft To gesetzt werden kann. Diese Vereinfachung 
erspart auch die komplexe Modellierung des Strahlungsflufies in der Atmosphare, die bereits mit 
einfachen Annahmen sehr umfangreich ist (78) . Bei Tagesfliigen fehlen ohnehin genaue Wetterda- 
ten, um den Einflufi der Wolkenbedeckung auf die Sonneneinstrahlung zu rekonstruieren. 



2 Die Nusselt-Zahl Nu ist proportional zum Verhaltnis der Warmeleitung einschlieBlich Konvektion und Turbulenz 
zur molekularen Warmeleitung. Es gilt Nu > 1. 
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2.5 Materialverhalten 



Bei ticfen Temperaturen verliert Gummi seine Flexibility und erstarrt in einem Glaszustand. Fiir 
einen aufsteigenden Ballon bedeutet dieser Phaseniibergang die Zerstorung der Hiille durch das 
sich weiterhin ausdehnende Fiillgas. Der Glasiibergang erstreckt sich iiber einen Temperaturbereich 
und lafit sich nur naherungsweise durch eine feste Temperatur beschreiben. Besonders fiir stark 
gedehnte Ballons mufi ein allmahlicher Ubergang vom Platzen durch Uberdehnung zum Platzen 
durch Verglasung angenommen werden |70] . 

Zur Modellierung dieses Effckts wird ein einfaches Modell verwendet: Unterschreitet die Ballon- 
temperatur eine Grenze Tpi atz , wird der Ballon als geplatzt betrachtet. Experiments von Kofoed- 
Hansen [55] zeigen deutlich, dafi bis — 20° C die elastischen Eigenschaften nicht beeintrachtigt 
werden. Burchette [TS] gibt die Glasiibergangstemperatur mit — 40° C an, wahrend Roberts [8~2"1 
einen Wert um — 60°C annimmt. 

Aus den in Kapitel [6] durchgefuhrten Simulationen, auf die an dieser Stelle vorgegriffen wird, 
ergibt sich die Notwendigkeit des Glasiibergangs bei einer Lufttemperatur von etwa — 55°C. Zu- 
sammen mit den in Kapitel 12.41 ausgefuhrten Uberlegungen ist nicht ausgeschlossen, daB der Bal- 
longummi eine andere Temperatur hat. In diesem Sinne ist die aus den Simulationen abgeleitete 
Temperatur hauptsachlich eine Information iiber die maximale Hohe, die ein Ballon unter opti- 
malen Umstanden erreichen kann. 



2.6 Aerodynamik 

Auf einen in z-Richtung aufsteigenden Gasballon der Masse M wirkt der Nettoauftrieb F a und 
die Reibungskraft Fr 

Mz = F a + Fji (2.23) 
wobei die Reibungskraft durch den Ansatz nach Newton beschrieben wird: 

F R = -^c w A q p L v z \v z \ (2.24) 

A q ist die Querschnittsflache in Bewegungsrichtung, c w der Widerstandsbeiwert und v z die verti- 
kale Geschwindigkeit des Ballons. Die Losung der Diffcrcntialgleichung 12.231 ist: 

z(t) = I In [cosh (tVoa/l)] (2.25) 
v z (t) = Woo tanh (tooo/i) (2.26) 

Die Endgeschwindigkeit Voa und die Langenskala I berechnen sich aus den Parametern des Ballons 
zu: 

_ 1 = (2.27) 
Mit einer Referenzluftdicht<|| ph,R e f kann diese Gleichung weiter vereinfacht werden: 




/PL,Rof / 2F a 

Woe = foo,oW , Vog,o ■= \ -. (2.28) 

V i°L y C w A q pL tRe { 

Ein Kartenballon erreicht innerhalb von Bruchteilen einer Sekunde die Endgeschwindigkeit Voq. 
Das Gesamtgewicht eines Ballons hangt von den verwendeten Karter@ und Materialien ab. 



3 PL,Ref = 1,204 kg/m 3 , siehe Anhang [XI 



4 Ballonpostkarten haben das iibliche Format 105 X 148 mm, sind aber mit 80 g/m 2 deutlich leichter als Stan- 
dardpostkarten. 
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Datum 


Durchmesser 


Hohe 


Zeit nach Start 




^z.thco 


Quelle 


30.08.1998 


6,6' 


24° 


45 s 


1,3 m/s 


1,32 m/s 


Photo 


20.09.1998 


21,0' 


31° 


15 s 


1,5 m/s 


1,34 m/s 


Photo 


21.07.1999 


15,6' 


24° 


25 s 


1,0 m/s 


1,51 m/s 


Photo 


05.08.1999 


9,0' 


38° 


40 s 


1,7 m/s 


1,45 m/s 


Photo 


14.01.2001 


ps 6" 


14° 


20 min 


1,8 m/s 


1,14 m/s 


Telcskop 



Tabelle 2.1: Auswertung der Steiggeschwindigkeit v z aus Ballonstarts. Die vorletzte Spalte enthalt 
theoretische Werte nach Gl. 12.271 Die Ballons sind Teil der in Kapitel [6T41 vorgestellten Versuchs- 
reihe. 

Eine kleine Auswahl bietet die folgende Tabelle: 



Postkarte 


3,8 g 


Ballonkarte 


1,6 g 


Geschenkband mit Plastikclip 


0,9 g 


Papierkarte in PE-Hiille 


2,0 g 



Fur den Widerstandsbeiwcrt c w findct man: 





Bemerkung 


0,44 


Litcraturwcrt Kugel [88] 


0,81 ±0,02 


Ballon alleine beim Aufsteigen 


0,93 ±0,09 


Ballon mit Karte 


0,96 


Pilotballon nach [591 



Wenn man die geringe Abweichung eincs Ballons von der Kugelform bedenkt, erscheint die deut- 
liche Abweichung des klassischen Literaturwerts fur eine Kugel zunachst erstaunlich. Der Unter- 
schied ist abcr nur zu einem geringen Teil auf die veranderte Geometrie zuriickzufuhren. Messungen 
des Stromungswiderstandes gehen von einem fixierten Priifkorper aus, was einem im Vergleich zum 
umgebenden Medium sehr dichten beweglichen Korper entspricht. Ein steigender Ballon hat aber 
offensichtlich eine geringere mittlere Dichte als Luft, so dafi die intensive Wechselwirkung mit der 
Wirbelschleppe den Luftwiderstand erhoht [5^1 [73] . 

Um diesen Befund unter realistischen Bedingungen weiter abzusichern, wurde der gefundene 
Cu,-Wert mit den Aufstiegsgeschwindigkeiten von angepeilten Ballons in Bodennahe in Tabelle HOI 
verglichen. Zu beachten ist, dafi die angegebenen Steiggeschwindigkeiten stark durch die Variabi- 
lity des Windes in Bodennahe schwanken (vgl. Abb. I6.2j) . 

Platzt der Ballon durch aufiere Einfliisse oder ein zu grofies Volumen, fallt die Karte zu Boden. 
Fiir eine einzelne Karte in einer PE-Hiille findet man Voo.o — 1, 0±0, 12 m/s. Die wirksame Flache 
A w = c w A q betragt in diesem Fall 326 cm 2 . 

2.7 Permeabilitat 

Die Gummihulle eines Ballons schliefit die Gasfullung nicht hermetisch ein, sondern erlaubt einen 
geringen Gasaustausch, der im Laufe der Zeit zum Schrumpfen des Ballons fuhrt. Der Gastransport 
durch die Gummihulle kann in zwei Prozesse zerlegt werden: Das Gas lost sich zunachst in der 
Membran (beschrieben durch das Henry-Gesetz) und diffundiert anschliefiend durch die Membran. 
Die Kombination beider Prozesse liefert die Stoffmenge pro Zeit N, mit der das Gas eine Membran 
der Flache A und der Dicke d durchdringt: 

A 

N = Q-^-Ap (2.29) 

Getrieben wird der Flufi durch eine Partialdruckdifferenz Ap. Aus diesem Grund ist der Gasverlust 
eines Luftballons deutlich kleiner als der eines Heliumballons, obwohl die Permeabilitaten selbst 
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5[10- H dm7h] 


melloc© 32 cm 


melloc© 13 cm 


Daten von [SB] 


Daten von 96 


N 2 


1,2 ±0,2 


1,2 ±0,7 




1,67 


o 2 


2, 2 ±0,9 


3, 3 ±1,7 




4,78 


Luft 


1,3±0,3 


1,6±0,7 


0,7 




H 2 


6, 2 ±0,5 


7, 4 ±1,0 


7,2 


10,71 


He 






3,8 


6,31 


co 2 




22 ±3 




27,53 



Tabelle 2.2: Permcabilitaten 5 des Ballongummis bei Zimmertempcratur und Normdruck. Die 
Volksbankballons sind mit den Daten der melloc© 32 cm Ballons konsistent. 



fur beide Gase vergleichbar sind. Die Permeabilitat Q laBt sich aus der Diffusivitat D und der 
Henry-Konstante fc# berechncn: 

Q = Dk H (2.30) 

Q ist temperaturabhangig. In Analogie zu Arrhenius-Gleichung gilt die Naherung 

Q(T) = Q exp(-T Q (l/T-l/T )) (2.31) 

wobei Tq die in eine Tempcratur umgerechnete Aktivierungsenergie der beteiligten Prozesse ist. 
Qo ist der bei der Temperatur To bestimmte Referenzwert. Fur die weiteren Betrachtungen wird 
Tq ss 1700 K (geschatzt nach [SB]) verwendet. Gleichung 12.291 lafit sich fur ein ideales Gas leicht 
mit Hilfe der Konzentration n schreiben: 

N = S-An (2.32) 
a 

Einsetzen der idealen Gasgleichung liefert die auf die Konzentrationsdifferenz An bezogene Grofie 

S := QRT (2.33) 
die weiterhin ebenfalls kurz als PermeabilitalU bezeichnet wird. Die Temperaturabhangigkeit ist: 

6(T) = S ^exp(-T Q (l/T-l/T )) (2.34) 
So := QoRTo (2.35) 

Die Modellierung der Temperaturabhangigkeit der Permeabilitat hat einen entscheidenden Einflufi 
auf den Gasverlust in grofier Hohe und damit die Flugweite. Tabelle [2^21 gibt einen Uberblick iiber 
die Permeabilitat von Ballongummi fiir verschicdene Gase. Die letzte Spalte enthalt Werte fiir 
vulkanisierten Naturgummi, die nur bedingt auf Ballongummi ubertragen werden konnen. 



2.8 Diffusion 

Mit den Permeabilitaten der beteiligten Gase kann der Volumenverlust eines Ballons als Funktion 
der Zeit bestimmt werden. In der Praxis sind zwei Gase - Luft und die Ballonfiillung - oft aus- 
reichend fiir eine genaue Beschreibung. Da Luft im wesentlichen ein Gemisch aus Stickstoff und 
Sauerstoff ist, bietet es sich aber an, noch ein weiteres Gas hinzuzufugen und einen Ballon mit drei 
verschiedenen Fiillgasen als Modell zu verwenden. Die Konzentration der drei Gase im Aufienraum 
sei unabhangig von der Zeit ni,n 2 und 713. Die Stoffmengen der Gase im Ballon werden durch 

5 Um Mehrdeutigkeiten auszuschlicficn wird in dicser Arbeit ausschliefilich <5 verwendet. 
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Ni, N 2 und N 3 gegeben. Jedes Gas hat eine individuelle Permeabilitat Si. Die Zeitentwicklung der 
Stoffmengen nach Gl. 12.321 ist darni: 

S 2 (n 2 -N 2 /V) (2.36) 
S 3 (n 3 -N 3 /V) J 

Die Dicke d der Gummihiille ergibt sich in sehr guter Naherung aus dem Gummivolumen der Hiille 
zu d — V<3 umm i/A. Alle weiteren Parameter wie der Innendruck Pi nnC n, der Aufiendruck PauBon und 
die Temperatur T werden als konstant angenommen. Zur Vereinfachung der Glcichung werden 
neue Variablen Vi — NiRT/p lnncn und r, = riiRT '/p au Sen eingefiihrt: 




A 2 1 



Si{p-n) -Sin 

-S 2 r 2 5 2 (p-r 2 ) -<i 2 r 2 | j V 2 | (2.37) 

-S 3 r 3 S 3 r 3 




VVcummi P 

i i 

p := El^>i (2.39) 

Paufien 

Zur weiteren Untersuchung wird die Matrix in Gl. 12.371 diagonalisiert. Die Eigenwerte 7* sind 
Losungen der Gleichung: 

Hieraus lafit sich direkt eine Ungleichung fur die Nullstellen ablesen: 

£3 > 73 > $2 > 72 > Si > 71 > (2.41) 
Die den Eigenwerten 74 entsprechenden Eigenvektoren sind 



1 

Vi = - 
P 



<5i — n 

S-2—li 

, S 3 r 3 
\ S 3 —ji 



(2.42) 



wobei jeder Eigenvektor ein Gasgemisch mit einem Einheitsvolumen beschreibt: 

3 

1 = ]T(^ (2-43) 

j=i 

Der Index j bezieht sich auf die einzelnen Komponenten des Vektors Sj. Da der Innendruck nur 
wenig grofier ist als der Aufiendruck, kann der kleinste Eigenwert 71 tiber eine Entwicklung nach 
p — 1 bestimmt werden: 

j " tt < 2M > 

l—l 

71 = S(p-l) + 0((p-l) 2 ) (2.45) 

6 kann als die mittlere Permeabilitat eines Gasgemisches angesehen werden. Fur den kleinsten 
Eigenwert 71 hat die durch v\ beschriebene Gasmischung praktisch die Gaszusammensetzung des 
Aufienraums. Dies entspricht einem luftgefullten Ballon, dessen Gasverlust ausschliefilich durch den 
Drucksprung der Ballonhiille getrieben wird. Die Diffcrcntialgleichungen fiir die Eigcnfunktionen 
Vi sind nun: 

= V+V 2 + V 3 (2.46) 
= - v f I 1 (2.47) 
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Einsetzen des Ansatzes 



Vi(t) = Vb,i exp(-tt(*)7<) (2.48) 
reduziert das Differentialgleichungssystem auf die Bestimmung der Hilfsfunktion u(t): 

* = ww^ < 2 ' 49 ' 

V(u)V Gumm i 

u{t) V(u) 

VGummi-rMdu (2.50) 

Gleichung 12.501 stellt implizit dar. Aus den Eigenfunktionen in Gl. 12.481 konnen dann die 
Funktionen Vi(t) mit Hilfe der Eigenvektoren in Gl. 12.421 linear kombiniert werden. 1st nur ein 
Vb,i 7^ 0, lafit sich die Losung explizit angeben: 

n ■■= —j2 — (2-51) 

Vi(t) = — o (2.52) 

(l + VOr,)) 3 

Zur Zeit t = Ti ist das Volumen etwa auf die Halfte abgefallen. Bei einem vernachlassigbaren 
Drucksprung (p = 1) hat das Diffusionsproblem einc Erhaltungsgrofie. Uber die Eigenvektoren Wi 
der transponierten Matrix 



(2.53) 



gewinnt man mit 71 = die explizite Darstellung: 

Y^Yl^l = k onst. (2.54) 

2—1 

Die obenstehenden Gleichungen konnen sofort auf hoherdimensionale Matrizen desselben Typs 
verallgemeinert werden. 
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Kapitel 3 



Wetter 



3.1 UV-Strahlung und Ozon 

Die Sonneneinstrahlung in Verbindung mit Oxidantien in der Luft hat einen starken EinfluB auf 
die Haltbarkeit eines Gummiballons. Die UV-Strahlung zerstort dabei den Gummi nicht nur direkt 
durch das Aufbrechen der Polymcrkcttcn, sondern vor allem durch die dadurch ermoglichten ra- 
dikalischen Reaktionen mit Ozon und Stickstoffoxiden. Einmal gestartet, laufen diese Reaktionen 
auch im Dunkcln wciter und fiihren schlieBlich zum Zerfall des Gummis bzw. Platzen des Ballons 

m 

Experimente von Burchette |15j bestatigen dieses Szenario, in dem durch Witterungseinfliisse 
und Sonnenlicht bereits im Laufe eines Tages der Grofiteil der Testballons zerstort wird. Die 
Zerfallsrate nimmt mit der Spannung des Gummis zu, so dafi besonders stark aufgeblasene Ballons 
zuerst platzen [68] . 

In milder Form macht sich die Oxidation des Gummis bereits nach einigen Stunden durch eine 
Mattierung der Ballonoberflache bemerkbar, die von einem erhohten Gasverlust begleitet wird 

Wahrend die Wichtigkeit der UV-Strahlung aufier Frage steht, ist die realistische Modellierung 
weitaus schwieriger. Ballonspannung, die Hysterese der Ballonhulle und die Volumenanderung 
beim Aufsticg wechselwirken mit der Oxidation des Gummis und den wechselnden Wetterein- 
fliissen. Ohne kontrolliert durchgefiihrte Experimente ist es kaum moglich, eine theoretische Be- 
schreibung zu gewinnen. Aus diesem Grund wird auf die direkte Modellierung dieses Prozesses 
zugunsten einer empirischen Untersuchung der Ballonflugdaten verzichtet. 

3.2 Regen und Luftfeuchtigkeit 

Regen ist allgegenwartig. Rund 50 % aller Tage in Deutschland gelten als Regentage, haben also 
einen Gesamtniederschlag von mehr als 0,1 mm in 24 Stunden 50]. Die Wichtigkeit auch geringer 
Niederschlage fur die Simulation eines Ballons laBt sich verdeutlichen indem man die Wirkung 
eines Regens mit Niederschlagsrate w auf den Auftrieb mit dem Verlust des Fullgases vergleicht. 
Die jeweiligen Auftriebsverluste F a durch diese beiden Prozesse sind: 

^a.Rogon ~ ~9wA c Qp W (3.1) 

A 2 

■Fa.Diff ~ -9PhS— (3.2) 

y Gummi 

A e s reprasentiert die effektive Sammelflache des Ballons fur Regen und pw ist die Dichte von 
fliissigem Wasser. Die weiteren Variablen entsprechen der Nomenklatur aus Kapitel |2~51 
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Aufloscn von F aj R egen < F ai Diff nach w ergibt die Bedingung fur die Dominanz des Regens beim 
Auftriebsverlust : 

p L A 2 5 

w - — l - ^ ( 3 - 3 ) 

PW A ctt "G ummi 

Einsetzen plausibler Werte ergibt eine Grofienordnung von 0,14 mm/Tag, also ein sehr leichter 
Nieselregen. Eine solche Wasseriibernahme konnte bereits bei einem Flug durch Wolken erreicht 
werden. Obwohl Ballongummi einen hydrophoben Charakter hat, kann sich genug Wasser in Form 
grofier Tropfen auf der Ballonhiille sammeln, um den Auftrieb eines Kartenballons vollstandig zu 
neutralisieren und ihn zu Boden driicken. Eine iiber den Ballonquerschnitt A q integrierte Nieder- 
schlagsmenge von 0,1 mm ist dafiir bereits ausreichend. Zufalligerweise entspricht diese Nieder- 
schlagsmenge dem Mindestniederschlag von 0,1 mm, um einen Tag als Regentag zu bezeichnen. 
Folglich bringt jeder ausgepragte Niederschlag einen Ballon mit grofier Wahrscheinlichkeit sofort 
zu Boden. 

Leider sind frei verfiigbare Daten zu grobmaschig, um die lokal sehr stark schwankende Nie- 
derschlagsverteilung ausreichend aufzulosen [21 [71] . Die grofite fiir die Forschung freie Wetterda- 
tensammlung in Europa ist das European Climate Assessment & Dataset [50] . Leider schwankt 
die Dichte der verfugbaren Stationen sehr stark und ist in vielen europalschen Landern zu gering, 
zudcm sind die Niederschlagsdaten nur integriert pro Tag verfugbar. Versuchsweise wurden diese 
Wetterdaten iiber Thiessen-Polygone interpoliert, da die geringe Stationendichte die Verwendung 
eines einfachen und robusten Verfahrens erfordert, und mit den Trajektorienrcchnungcn (siehc 
Kap. [5]) kombiniert. Eine aussagekraitige Auswertung scheiterte aber an der geringen Auflosung 
der Daten, so dafi dieser Ansatz nicht weiter verfolgt wurde. 



3.3 Atmosphare 



Die Atmosphare gliedert sich in mehrere verschiedene dynamische Schichten. Direkt an den Boden 
grenzt die laminare Unterschicht mit einer Starke von wenigen Millimeter, die durch die Viskositat 
der Luft dominiert wird. Daran schliefit sich die Prandtl-Schicht an, ein turbulenter Luftstrom der 
in den unteren Bereichen von der Oberflachenstruktur bestimmt wird. Sie ist 20 bis 60 Meter 
dick. Die Eckman-Schicht bildet schliefilich den Ubergang zu dem durch Druckgradienten und 
Corioliskraft gepragten geostrophen Wind. 

Das bodennahe Windprofil ist fiir den Start eines Kartenballons von besonderem Interesse. 
Eine oft verwendete Beschreibung ist die Monin-Obukhov-Ahnlichkcitsthcorie. Sie basiert auf ei- 
ner dimensionslosen Formulierung der relevanten hydrodynamischen Gleichungen und einer ge- 
schickt gewahltcn Schliefiungsbedingung. Dieser Ansatz reduziert die Diffcrentialgleichung fiir das 
hohcnabhangige Geschwindigkcitsprofil v(z) auf 

f % = •«(«/ £ ) <") 

mit der von-Karman-Konstante K, der Reibungsgeschwindigkeit v» und der Obukhov-Lange L. 
Eine kompakte Darstcllung mit weiterfiihrender Literatur ist in |31j enthalten. Die verbleibende 
universelle Skalenfunktion $m mufi experimentell bestimmt werden [16] . Fiir neutrale Schichtun- 
gen ist L — > oo und die Skalenfunktion reduziert sich auf die Konstante = 1. In diesem Fall 
lafit sich Gl. 13.41 einfach integrieren und man erhalt das klassische logarithmische Windprofil: 

v(z) = ^Hz/z ) (3.5) 

Die Integrationskonstante zq ist die sogenannte Rauhigkeitslange, die von der Gelandestruktur 
abhangt und ebenfalls experimentell ermittelt werden mufi. 

Ein erster Schritt, um die Atmosphare in grofien Hohen zu beschreiben, sind gemittelte Werte 
von Wind, Druck und Temperatur als Funktion der Hohe. Als Grundlage werden die Daten aus 
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Tabelle lE.il fur die mittlere Troposphare iiber Deutschland verwendet, urn ein einfaches und rea- 
listisches Modell zur Hand zu haben. Diese Daten lassen sich gut durch folgende Funktionen in 
Abhangigkeit der Hohe z darstellen: 

p[Pa] = 101325exp^-^— ) (3.6) 

v[m/sl = 2,58— ^— + 5,13 (3.7) 
1 km 

T[°C] = To- 4, 223^ - 0, 187 " (3.8) 

p[kg/m 3 ] a i )2 6expf-g-^— ) (3.9) 

Die Temperatur To am Boden ist im Mittel 9°C, kann aber verandert werden urn verschiedene 
Jahreszeiten zu reprasentieren. Fiir einfache Studien des Flugverhaltens sind diese Gleichungen 
ausreichend und werden im weiteren als Referenz verwendet. Die Hohenabhangigkeit der Dichte 
kann aus dem Druck- und Temperaturverlauf berechnet werden und ist hier nur fiir analytische 
Abschatzungen mit angegeben. 

Die genaue Berechnung von Trajektorien erfordert zwangslaufig bessere zeit- und ortsauf- 
geloste Daten fiir den betrachteten Zeitraum und Startorte. Als Datenquelle kommen praktisch 
nur Radiosondenaufstiege in Frage, die die gewiinschten Werte punktuell als Funktion der Hohe 
liefern |67j . In dieser Arbeit werden Daten von [79 verwendet, die kostenfrei fiir Forschung und 
Lchre zur Verfiigung stehen. Die im HTML-Format vorliegenden Daten wurden konvertiert und 
einem einfachen Plausibilitatstcst unterworfen, um Datenfehler zu erkennen. Die Interpolation in 
der Hohe erfolgt mit einer Spline- Interpolation dritter Ordnung. Datensatze mit zu wenigen Punk- 
ten wurden verworfen, um Interpolationsproblcme zu vermeiden. Die raumliche Interpolation der 
Daten Xi an den Stiitzstellen Xi erfolgt mit einer einfachen inversen Distanzgewichtung |47j 

X (x) = Ik^i (3.10) 

m = J^W (3 ' n) 

wobei X ein Platzhaltcr fiir die zu intcrpolierende Grofie am Ort x ist. Abweichend von dem oft 
verwendeten Exponenten 2 in der Gewichtungsfunktion wurde der Exponent 4 gewahlt, der ein 
glatteres Geschwindigkeitsfeld liefert und in Vergleichen mit den Ballonflugdaten geringere Fehler 
aufzeigt. Besonders kiinstliche Maxima der Windgeschwindigkeit an den Sttitzpunkten („Bullseye- 
Effekt") werden reduziert. Der Ansatz [5.101 gewahrleistet. dafi die interpolierte Grofie nicht das 
Maximum der diskreten Werte iibersteigen kann. Dies gibt dem Verfahren eine grofie Robustheit, 
bedeutet aber auch eine systematische Unterschatzung von z. B. interpolierten Windgeschwindig- 
keiten. 

Neben der inversen Distanzgewichtung stehen eine Vielzahl von Interpolationsalgorithmen zur 
Verfiigung (siehe die Ubersicht in [55]), die zum Beispiel statistische Korrelationen oder ausgefeil- 
tere Basisfunktionen verwenden [381 1103] und dadurch eine hohere Genauigkeit erreichen konnen. 
Da die Interpolationsqualitat von Radiosondcndaten bereits durch die Dichte der Messpunkte 
limitiert ist [93] und komplexe Modelle durch die fehlende Verfiigbarkeit grofier Rechenleistung 
von vornherein ausgeschlossen waren, wurde die inverse Distanzgewichtung als bester Kompromiss 
gewahlt. Die zeitliche Interpolation der Daten erfolgt mit einer einfachen linearen Interpolation, 
um storende numerische Artefakte bei dynamischen Wettersituationen zu vermeiden. 

Eine anschauliche Darstellung der jahreszeitlichen Effekte gibt die Datenstatistik der Stati- 
on Meiningen in Abbildung 13.11 Das Temperaturprofil folgt der Temperatur am Boden, variiert 
aber sonst nur wenig. Der Gang der Windgeschwindigkeit mit der Jahreszeit zeigt eine halbjahr- 
liche Periode mit einem starken Maximum im Winter und einem weiteren Maximum im Sommer, 
die auf Strahlstrome in rund 10 km Hohe zuruckzufuhrcn sind. Die noch deutlich erkennbaren 
Schwankungcn des gcmitteltcn Windcs dcutcn auf einen zu kurzen Mittclungszcitraum hin. 
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Windgeschwindigkeit 10548 Meiningen 




Monat 




1 2 3 4 5 6 7 8 9 10 11 12 



Monat 



Abbildung 3.1: Uber die Jahre 2000-2008 gemittelter Wind- und Temperaturverlauf der Radioson- 
denstation ,,10548 Meiningen", geglattet mit einem GauB-Kern der Breite 7 Tage. h max markiert 
die — 55°C-Isotherme. 
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Die Fluktuation der Windgeschwindigkeit um dicscn Mittelwert wird gut durch eine zweipa- 
rametrige Weibull-Dichtcfunktion wiedergegeben [72] : 

f(v) = A^"^-^)*} (3.12) 

Der Formparameter fc ist vom Ort, der Hohe und der Jahreszeit abhangig. Der jahresgemittelte 
Wert liegt im Bereich fc = 1, 8 ... 2, 6, so dafi fc = 2 eine realistische Wahl ist. Fiir fc = 2 entspricht 
die Weibull-Verteilung der Rayleigh-Verteilung. Mittelwert E und Varianz Var sind: 

E(v) = u r(l/fc + l) (3.13) 

= vq^- fiir fc = 2 (3.14) 

Var(v) = vl (r(2/fc + l) -r(l/fc + l) 2 ) (3.15) 
= «g(l-7r/4) fiir fc = 2 (3.16) 

r ist die bekannte Gammafunktion. 



3.4 Turbulenz 

Neben den groBskaligen Luftbewegungen, die als Wind wahrgenomnien werden, ist jede Luft- 
stromung auch von Turbulenz bis hin zu den kleinsten Skalen gepragt. Fiir die Simulation eines 
Ballons ist es sinnvoll, zwei verschiedene Anteile der Turbulenz getrennt zu betrachten: Vertikale 
Luftbewegungen, die das Aufsteigen des Ballons beeinflussen, sowie Fluktuationen der horizon- 
talen Driftbewegung. Senkrechte Luftbewegungen beeinflussen das Aufsteigen des Ballons und 
entscheiden somit, wann und ob der Ballon seine Platzgrenze erreicht. Dariiber hinaus kann ei- 
ne veranderte Verweildauer in den verschiedenen Luftschichten bei starken Windscherungen auch 
die Flugrichtung beeinflussen. Das Ausmafi vertikaler Luftbewegungen hangt wesentlich von der 
Stabilitat der Luftschichtung ab: Nur bei labilen Wetterlagen und starker Konvektion, in starken 
Windscherungen oder in der Nahe von orographischen Hindernissen ist sie bemerkbaiQ. Im Mit- 
tel ist die vertikale Luftbewegung recht gering; Messungen ergeben typische Geschwindigkeit im 
Bereich weniger Meter pro Sekunde [75 . Verkehrsflugzeuge sind nur wenige Prozent der Flugzeit 
Turbulenzen ausgesetzt|29). was ebenfalls fiir eine ruhige Atmosphare spricht. Gegeniiber die- 
ser typischen Situation, in der die vertikale Turbulenz auf den Flug eines Ballons vernachlassigt 
werden kann, ist der Einflufi bei starker Konvektion unbestreitbar. Da diese Situation oft von 
Niederschlagen und Gewittern begleitet wird, ist ein abruptes Ende des Ballonflugs aber wahr- 
scheinlicher als eine systematische Bccinflussung der Trajektorie. Aus diesem Grund wird dieser 
Effckt vernachlassigt und spater iiber die empirischen Daten erschlossen. 

Der horizontale Anteil der Turbulenz beschrankt nur die erreichbare Genauigkeit der Trajek- 
torie. Eine erste Quantifizierung der Turbulenz bietet die Korrelation der Geschwindigkeiten v an 
zwei Orten x\ und x 2 zu einer festen Zeit t: 

C v (xi,x 2 ) ■= (v(xi) ■ v(x 2 )) (3-17) 

Das Mittel ist hier iiber verschiedene Realisierungen der Stromung zu fiihren. Fiir die Betrachtung 
von Trajektorienfehlern ist es sinnvoller, die Korrelation der Geschwindigkeitsdifferenz Av zu 
betrachten (Strukturfunktion zweiter Ordnung): 

Av := v 2 - vi (3.18) 
(Av-Av) = C v {x 1 ,x 1 ) + C v {x2,x 2 )-2C v {x 1 ,x 2 ) (3.19) 

1 Die Stratosphare zeigt durch die ihr eigene Tempcraturinversion eine besonders groBe Stabilitat. Ein Karten- 
ballon verweilt aber praktisch nur in der Troposphare. 
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Bei raumlicher Homogenitat der Stromung entfallt die absolute Ortsabhangigkcit: 

(Aw- Av) = 2(v 2 ) - 2C v (x 2 - xi) (3.20) 

Eine alternative Charakterisierung ist die spektralen Verteilung der kinetischen Energie, die aus 
dem Mittel der fouriertransformierten Geschwindigkeiten v(k) gewonnen wird: 

E(k) := ^(|^)| 2 > (3.21) 

Die Geschwindigkeitskorrelation steht iiber dem Faltungstheorem in Zusammenhang mit der spek- 
tralen Verteilung der kinetischen Energie 

E(k) = \c v (k) (3.22) 

wobei C v (k) die fouriertransformierte Korrelation aus Gl. 13.171 ist. Bei raumlicher Isotropic kann 
das Spektrum weiter zu einer eindimensionalen Funktion reduziert werden. Fur ein zweidimensio- 
nales Problem (z. B. stabil geschichtete Atmosphare) erhalt man die Projektion 

E{k) = J E{k')8{k-\k'\)dk' 2 (3.23) 

wobei hier 5 die Dirac-Funktion ist. Die Strukturfunktion und das Energiespektrum der Atmo- 
sphare sind zeit- und ortsabhangig. Globale Mittelwerte sind aber ein robustes Mafi fur die Verhalt- 
nisse in der Atmosphare [85] ■ Als Basis fur die weitere Betrachtung wird deshalb das global ge- 
mittelte Ergebnis von [ST] fur eine Hohe von etwa 10 km angegeben: 

(Av- Av) [m 2 /s 2 ] = 7,6 x 10~ 3 r 2/3 + 8, 9 x 10 _ V - 0, 59 x HTVlnr (3.24) 
E(k) [m 3 /s 2 ] = 9, 1 x 10~ 4 £;- 5/3 + 3 x lO" 10 /^ 3 (3.25) 
r = \2 2 -2i\ (3.26) 

r und k sind jeweils in SI-Einhcitcn cinzusetzen. 

Die in Gl. I3.24l beschriebenen Geschwindigkeitsfluktuationen treiben zwei benachbarte Ballons 
mit der Zeit auseinander. Diescr Effekt ist dabei fundamentaler Natur: Selbst wenn die Startor- 
te der Ballons beliebig nahe beieinander sind, entfernen sich die Ballons durch diesen diffusiven 
Prozcss. Diescr Aspekt deterministischen Chaos limitiert die Genauigkeit einer Trajektorienbe- 
rechnung unabhangig von der Datenqualitat und dem verwendeten Algorithmus. Fur die mittlere 
Entfernung r gleichzeitig gestarteter Ballons als Funktion der Zeit t gilt das klassische Ergebnis 
von Richardson [5D] 

(r 2 it)) = ?|(fe) 8 (3.27) 
= 0,6 2 ^; 3 (3.28) 

mit einem Wert des Parameters e = 0,0187 m 2 / 3 /s. Die Giiltigkeit dieses Gesetzes hangt von 
der Luftschicht und der absoluten Grofie der Ballonverteilung ab. Auf globalen Skalen geht das 
Richardson- Gesetz in eine normale Diffusion iiber |54) . 

Die mittlere Entfernung in Gl. 13.271 gibt nur eine grobe Vorstellung fiber die Verteilung der 
Ballons im Raum. Eine genauere Charakterisierung erlaubt eine Dichtcfunktion f r (r,t), welche 
die Verteilung der Ballons als Funktion der Zeit t und dem Abstand r zum Mittelpunkt der 
Ballonverteilung angibt. Uber die genaue Form dieser Dichtefunktion herrscht weniger Gewifihcit 
(vergleiche [3 HI [TUl [H]). Allgemein sind die Losungen aber von der Form f r oc exp(— r a ) [TTj . 
wobei es verschiedene Ansichten iiber die richtige Wahl des Exponentcn a gibt. Richardson gibt 
die eindimensionale Dichtcfunktion 

9 / 9lrl 2 / 3 \ 

''<'■«> " v^ exp (-J^ (3 ' 29) 
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in seiner Arbeit an. Bei einer maximalen Ballonflugzeit von 20 Stunden liefert Gl. l3.28l eine unver- 
meidbare mittlere Streuung in den Trajektorien gleichzeitig gestarteter Ballons von rund 60 km. 
Auch wenn dieser Wert eine kilometergenaue Berechnung von Ballontrajektorien ausschliefit, ist er 
deutlich kleiner als die Strecke von mehreren 100 Kilometern, die ein Kartenballon typischerweise 
in 20 Stunden zuriickgelegt. Eine Vertiefung dieser Fehlerbetrachtung erfolgt in Kapitel HTSl 
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Kapitel 4 



Boden 



Die Fundwahrscheinlichkeit einer Ballonkarte hangt davon ab, in welcher Umgebung der Ballon zu 
Bodcn gcht. Um einen Einblick in ein typisches Landegebiet zu erhalten, wurde mit dem Programm 
GoogleEartff] ein Experiment durchgefuhrt. Aus einem Quadrat^, das durch den 49ten und 50ten 
Breitcngrad sowie den 8,5°(E) und 9,5°(E) Langenkreis begrenzt wird, wurden 70 Landeplatze 
per Zufall ausgewahlt. Die Zusammcnsctzung der gcwahlten Flache aus Wald, Landwirtschaft 
und Siedlungen entspricht in etwa dem zu erwartenden Zielgebiet. 

Fur jeden fiktiven Fundort wurde der Abstand zum nachsten Weg (d. h. Strafie/ Feldweg/ 
Waldweg/ . . . ) und zum Rand der nachsten Siedlung sowie der Typ (Stadt/Wald/Feld) der Flache 
bestimmt. Die Anteile der verschiedenen Flachentypen sind in Tabelle B~T1 zusammengefafit. Die 
aus dieser kleinen Stichprobe ermittelten Daten passen gut zu den umfangreicheren Erhebungen 
aus [44], die ebenfalls in der Tabelle enthalten sind. 

Fur die Verteilung des Abstandes r zur nachsten Strafie bzw. Stadt wird eine kurze theoreti- 
sche Betrachtung durchgefuhrt. Abbildung 14.11 gibt eine schematische Zusammenfassung der im 
Folgenden verwendeten Groficn. Bchndct sich eine einzelne Struktur i mit der Flache As,i in einem 
Gebiet mit der Flache Ao, ist die Wahrscheinlichkeit, sich im Abstand r oder weiter entfernt zu 
bcfindcn 

P(>r) = (l-^-) (4-1) 



A ~ A s ,i 

wobei die Nachbarflache Ai der Bereich aller Punkte ist, die naher als r an der betrachteten Struk- 
tur sind. Durch zufalliges Hinzufugen weiterer (nicht notwendigerweise identischer) Strukturen bei 
glcichzcitiger Erweiterung der Gesamtflache erhalt man 



i=l 



Der Grenziibergang n — > oo fuhrt schliefilich auf die Poisson- Verteilung (vgl. z. B. [18] ) 

p(>r) = exp(-SA(r)) (4.3) 

n 

AM := lim -V^fr) (4.4) 



n— >oo TL 

i=l 



mit der neu eingefiihrten Anzahldichte der Strukturen S = 1/Aq. Die Abstandsabhangigkeit der 
Nachbarflache steht in direktem Zusammenhang mit der Hausdorff-Dimension Dh des Randes der 
betrachteten Struktur: 



A(r) cx r 2 -° H (4.5) 



^iehe http://earth.google.de/index.html 20.11.2006 

2 Das Quadrat wird in etwa durch die Stadte Karlsruhe — Frankfurt - Lohr am Main - Sulzbach an der Murr 
begrenzt. 
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Abbildung 4.1: Darstellung des kiirzesten Abstandes r vom Fundort zur nachsten Strafie und zum 
nachsten Ort. Die in Gl. 14.11 verwendeten Flachen sind fiir einen Fundort eingezeichnet. 



Umgekehrt kann der Langenskala r eine fraktale Dimension zugeordnet werden: 

Nach diesen Vorarbeiten werden nun die Abstandsverteilungen zur nachsten Strafie und zum 
nachsten Stadtrand empirisch aus dem Experiment ermittelt: 

PWcs (>r) « exp(-^-(l + y^Z) ) ,< 600 m (4.7) 

Pstadt(> r) cxp ^- Q ^— ^ J r < 2 km (4.8) 

Die Giiltigkeit dieser Verteilungen ist durch den geringen Umfang der Stichprobe auf die ange- 
gebenen Bereiche beschrankt. Wahrend die funktionale Form dieser Verteilungen rein empirisch 
bestimmt und daher ohne besondere Bedeutung ist, lafit sich der geometrische Gehalt gut mit 
Hilfc von Gl. 14.61 crschlicfien. Strafien erscheinen auf kleinen Skalen eindimensional, gehen aber 
auf grofieren Skalen in ein Netzwerk mit Dimension 1,5 iiber. Stadte erscheinen als aufgeloste 
Punktmenge mit Dimension 0,5, was gut den landlichen Charakter des bctrachtctcn Gcbictcs ab- 
bildet. Dieser Bcfund pafit zu weitergehenden Untersuchungen zu den fraktalen Eigenschaften von 
Siedlungsstrukturen [221 133] . die Strafien Dimensionen um 1,5 und landlichen Siedlungsstrukturen 
Dimensionen kleiner 1,0 zuordnen. 

Wenn man die regelmafiig von Menschen aufgesuchte Flache (d.h. die Flache, auf der ein Kar- 
tenballon hauptsachlich gefunden wird) mit den Angaben zu Siedlung und Verkehr in Tabelle 14.11 
gleichsetzt und problematische Fundorte wie Dacher und Baume vernachlafiigt, ergibt sich eine 
obere Grenze der Fundquote von etwas mehr als 10 %. In diinn besiedelten Gebieten, z. B. in 
Gebirgen, weiten Teilen der USA oder Osteuropa ist dieser Wert sicher noch nach unten zu korri- 
gieren |91) . Dariiber hinaus konnen Lander- bzw. Sprachgrenzen den Anted der zuriickgeschickten 
Karten merklich reduzieren (pers. Mitteilung M. Kollefrath). 
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[%} 


B.-Wurt. (2001) 


D(2001) 


GooglcEarth 


Landwirtschaft 


46,8 


53,5 


46 ±6 


Wald 


38,0 


29,5 


44 ±6 


Wasser/Sonstiges 


2,0 


4,3 




Siedlung u. Vcrkehr 


13,2 


12,3 


10 ±4 


davon: Gebaudeflache 


53,2 






Erholungsflache 


5,0 






Verkehrsflachc 


40,2 







Tabelle 4.1: Flachennutzung in Baden- Wiirttemberg und Gesamtdeutschland nach [33]. Die letzte 
Spaltc cnthalt das GoogleEarth-Experiment. Wald bezieht sich hier auf baumbestandene Flachen. 
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Kapitel 5 

Simulationen 



5.1 Gleichungen 

Die in den Kapiteln[5]-[3] vorgestellten Grundlagen bilden einen vollstandigen Satz von Gleichun- 
gen, um die Trajektorie eines gasgefiillten Ballons zu berechnen. Die Position des Ballons wird 
durch die geographische Lange A, die Breite /? und die Hohe iiber Normalnull z beschrieben. Die 
Windgeschwindigkeit ist durch v u (ostwarts) und v v (nordwarts) gegeben, vertikale Luftbewegun- 
gen werden vernachlafiigt. Die Massenbilanz und der Auftrieb sind: 

M = M + m G N G + m L N L (5.1) 
F a = Vg Ph (\,p,z)-Mg (5.2) 

Die Ballonmasse M ergibt sich aus dem Hiillen- und Kartengewicht Mo, der Stoffrnenge N G des 
Fiillgascs „G" und der eindifTundicrtcn Luftmenge Nl. m G und tul sind jeweils die Molmasse des 
Fiillgases und der eindiffundierten Luft. Die Trajektoriengleichung x — v(x, t) in den gewahltcn 
Koordinaten ist: 

P = (5 .3) 

^ = V U {\P,Z) 

COs(/3)(i? E rdc + z) 

^ ^Iw^) Ballon intakt. 

sonst - 

Das Ballonvolumen V und die Querschnittsflache A q werden mit der Zustandsgleichung des Bal- 
longases ermittelt: 

_ (N L +N G )RT L (\,p,z) 

PL (\,f3,z) + Ap { -> 
A = 5,0V 2/3 (5.7) 

A q = ^ (5.8) 

Pl und Tl sind das zeit- und ortsabhangige Druck- und Tempcraturprofil der Atmosphare (siehe 
Kap. 13.31) . Zuletzt wird die Diffusion des Fiillgases bestimmt: 

A 2 N r 
N G = - v d G (\,/3,z)-^ (5.9) 

'Gummi * 

Nl = --^—S L (X,/3,z)(^-n L (X,l3,z)) (5.10) 
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Ap 


5300 


Pa 


'max 


30 


1 


M 


7,1 


g 


^Gummi 


5,1 


g 




0,93 






1,0 


m/s 




-55 


°C 



Tabelle 5.1: Referenzparameter zur Beschreibung eines Kartcnballons. 

Der Ballon gilt als geplatzt, sobald das Volumen grofier als das Maximalvolumen V max ist oder die 
Temperatur unter die Glasiibergangstemperatur Tpi atz des Gummis fallt. 

5.2 Theorie 

Zunachst wird eine analytische Betrachtung der Gleichungen durchgefiihrt, urn die grundlegen- 
den Effekte zu verstehen. Zu diesem Zweck wird das rein hohenabhangige Atmospharenmodell in 
Gl. I3.61f3~8l und ein kartesisches Koordinatensystem fur die Ballonposition verwendet. Die Glei- 
chungen [52H53] andern sich dadurch auf: 

x = v(z) (5.11) 
V = (5.12) 

Fiir die folgenden Studien wird ein typischer Ballon betrachtet, dessen Parameter in Tabelle 15.11 
zusammcngestellt sind. Dieser Paramctcrsatz wird auch in Kapitel |5] verwendet, wenn keine de- 
taillierten Informationen iiber die gestarteten Ballons vorlagen. Durch das mogliche Platzen des 
Ballons wahrend des Fluges konnen drei verschiedene Trajektorientypen unterschieden werden: 

I Der Ballon bleibt wahrend des gesamten Fluges intakt. 

II Der Ballon erreicht die Maximalhohe und platzt durch den Glasiibergang des Ballongummis. 

Ill Der Ballon bleibt unter der Maximalhohe, platzt aber durch die Uberdehnung der Hiille. 

Fall U ermoglicht offensichtlich die grofiten Flugweiten, so dafi dieses Szenario am detailliertesten 
diskutiert wird. Anschliefiend werden die beiden letzten Szenarien angesprochen. 

Fall[D 

Um ein ubersichtlicheres Gleichungssystem zu erhalten, wird die Diffusion der Luft sowie der 
Drucksprung des Ballons vernachlaBigt (5l — 0, Ap = 0). Das Eindiffundieren von Luft vergroBert 
den Ballon und erhoht dadurch indirekt den Verlust des Fiillgases. Dariiber hinaus verringert sich 
fiir Ap > der Auftrieb, was aber nur in grofien Hohen einen merklichen Beitrag liefert. Insge- 
samt ist der Einflufi auf die Ballontrajektorie aber zu vernachlassigen. Der endliche Drucksprung 
des Ballons sorgt fiir einen zusatzlichen Verlust des Auftriebs beim Aufsteigen, der aber „gespei- 
chert" und beim Abstieg wieder freigesetzt wird. Durch die starkere Kompression in grofien Hohcn 
kann der Ballon am Boden auch viel starker gefiillt werden, ohne zu platzen. Dies fiihrt zu einer 
deutlichen Asymmetrie der Trajektorie und ermoglicht grofierc Flugweiten. Trotz des dcutlichen 
Einflusses des Drucksprungs auf die Trajektorie, der in Kapitel I5~3l genauer untersucht wird, gibt 
eine vereinfachte Betrachtung wertvolle Einsichten. 
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Unter Verwendung dieser Naherungen erhalt man aus Gl. 15. II 15.21 und Gl. I5.6H5.121 

F a = {N G (m L - m G ) - M )g (5.13) 

z = sign(F a ),/ 2 f a \ \ ( 5 ' 14 ) 

y C w AqP L {Z) 

x = v[z) (5.15) 
Ng = -tt^—5 g {z)^ (5.16) 

Zunachst werden die Gleichungen zusammengefaBt und dimensionslos formuliert. Aus der Leer- 
masse Mq, der Luftdichte /?l und der Permeabilitat Sq am Boden konnen alle notigen SkalengroBen 
gewonnen werden: 

M o 

V := t := — —2 — (5.17) 

/ := -T7— v z , :=\ (5.18) 

M a g y c w A p Lfi 

Vq ist das Volumen eines Ballons mit F a — 0, tq die entsprechende Gasverlustzeitskala, / der 
Auftrieb in Einheiten des Gesamtgewichts und v z ^ die Steiggeschwindigkeit eines Ballons mit 
Auftrieb / = 1. Die Luftdichte p und Permeabilitat 8 wird in Einheiten der Werte am Boden 
(z — 0) ausgedriickt: 

p ■= my 6 - = W) (5 9) 



Die dimensionslosen Gleichungen sind: 

/ = 



1(1 + /) 4/3 ^ (5.20) 

TO \/ph 

(l + /)V3^/« 



sign(/H,o- Vlf [ 1/fi (5.21) 



Dicsc konnen zu 



sign(/) 7i 7 «K?a <y ( 5 - 22 ) 



J/e L * wy (l + /)»/3 

:= «z,oTo (5.23) 

zusammengefaBt werden. Die linke Seite enthalt nur von 2; abhangige Grofien und kann integriert 
werden: 

* G dZ ' (5.24) 



p^ L Q 

In erster Naherung kann z sa z/Lq gesetzt werden. Zunachst wird / <C 1 untersucht. Die Integra- 
tion der Differentialgleichung ab z = liefert: 

S(f) = l(fT-\f?' 2 )+0{fl*) (5.25) 

Die maximale Flugweite wird erreicht, wenn der Ballon knapp unter der Maximalhohe h max bleibt 
und intakt zu Boden sinkt. Damit folgt fur den optimalen Auftrieb: 

/3~ \ 2/3 

F a ,o P t = M g (5.26) 
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Wasserstoff, V max = 30 1 





Abbildung 5.1: Ballonfliige in Abhangigkeit von der Leermasse Mo und dem Auftrieb F a . Die 
Konturlinien beschreiben die jeweils erreichte Scheitelhohe h* . Zum Vergleich sind die Ballons aus 
Tabelle 16.21 mit eingetragen. 
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Die Flugweite ergibt sich aus der Zeitintegration des Windprofils in Gl. 13.71 



d = / v(z(t))dt (5.27) 
Jo 

« 2r Q f v(3/5h nmx ) (5.28) 



Fiir den optimalen Startauftrieb nach Gl. 15.261 erhalt man: 

d max « ,(3/ 5 ^) ( %^ 2/9 ( /^^ f _^_y /9 (5 29) 

V M) / y ff<5G,o \m L -m G J 

rfmax ist die maximale Flugweite, die mit einem gegebenen Kartenballon erreicht werden kann. Bei 
grofieren Werten der Leermasse Mq ist auch ein grofierer Ballon notig, so dafi der Ballon schliefilich 
nicht mehr die maximale Hohe h max erreichen kann und vorher platzt (Ubergang zu Fall IHip . In 
diesem Fall muB erst die Scheitelhohe /i* ermittelt werden: 

V max p L (h*) = M(h*) (5.30) 

Mit dieser kann dann der optimale Auftrieb entsprechend berechnet werden: 

3 X 2/3 



Fallini 



F a , op t = M g [ -K ) (5.31) 



Wird der Auftrieb geringfiigig gesteigert, platzt der Ballon durch die Kalte und fallt zu Boden. Die 
nun erreichte Entfernung dtmrst die maximal Flugweite fiir einen Ballon, der vor dem Erreichen 
des Scheitelpunktes platzt und zu Boden fallt. Dieser Fall tritt ein fiir: 

< F a (h max )/g < V max p L (h max ) - M(h max ) (5.32) 

h* = hmax (5.33) 

In diesem Regime ist die Scheitelhohe gleich der Maximalhohe. Die Ungleichung 15.321 kann fiir 
Kartenballons mit grofier Leermasse nicht erfiillt werden. In diesem Fall ist das Erreichen der 
Maximalhohe unabhangig von der Ballonfiillung ausgeschlossen. 



Fall HI] 

Ab einem bestimmten Auftrieb platzt der Ballon bereits durch die iiberdehnte Hiille, bevor die 
Kalte in grofien Hohen den Ballon zerstort. Dies definiert eine weitere Flugweite dburst2- Der Ballon 
platzt in der maximalen Hohe /i max (oder friiher) wenn 

F a (h max )/g > V max ph(h max ) -M(/i max ) (5.34) 

gilt. Fiir den dimensionslosen Auftrieb in dieser Hohe bedeutet dies: 

f(h max ) > V max ph{h max ) / M - M(h max )/M (5.35) 

Im Bereich grofier Auftriebe (Ubergang zu FalllHl) wird Gl. I5.25l ungenau und wird durch folgende 
Naherung ersetzt (mit einem Fehler von 10 % fiir 0, 2 < / < 3): 

2 « l(fo-f) (5.36) 
Nach dem Einsetzen dieser Naherung erhalt man fiir den Startauftrieb am Boden die Ungleichung: 

F a ,o > M max g - (1 - 3h max )M g (5.37) 

M max RJ V max p L (/lmax) (5.38) 
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Abbildung 5.2: Flugweite als Funktion des Auftriebs fur zwei Gasfiillungen. 

Dies ist eine notwendige Bedingung fiir das Platzen des Ballons vor dem Erreichen der Maxi- 
malhohe. Nur Ballons, die leichter als M max sind, konnen die Maximalhohe intakt erreichen. Der 
genaue Wert hangt vom Fiillgas ab und liegt im Bereich von 10 - 11 g (vgl. Abb. I5.1|) . Die 
Scheitelhohe in diesem Regime ist die Losung der Gleichung: 

F a {K)/g = V maxPL (h*)-M(h*) (5.39) 

5.3 Numerische Integration 

Zur genaueren Untersuchung wurden die Modellgleichungen mit einem Runge-Kutta Verfahrcn 
zwcitcr Ordnung numerisch integriert [1]. Als Zeitschritt wurden 50 Sekunden gewahlt. Durch 
diesc Wahl konnen Integrationsfehler vernaclilafiigt wcrdcn und die crzielte Hohcnaunosung von 
rund 100 Metern entspricht gut der vertikalen Auflosung der interpolierten Wetterdaten. 

Die in der Praxis wichtigsten (und am leichtesten beeinflufibaren) Parameter eines Kartenbal- 
lons sind die Leermasse, der Startauftrieb (bzw. die Ballonfiillung) und die Wahl des Fiillgases. 
Abbildung 15.11 zeigt die Integration des Musterballons aus Tabelle 15.11 jeweils mit Helium und 
Wasserstoff als Gasfullung und verschiedenen Leermassen und Startauftriebskraften. Zur deut- 
licheren Darstellung sind auch Leermassen kleiner 5,1 g enthalten, obwohl diese einer Nutzlast 
mit negativer Masse entsprechen. Zur besseren Vergleichbarkeit wird die Auftriebskraft mit der 
Eixlbeschleunigungfj] in eine aquivalente Masse umgerechnet dargestellt. Die bereits beschriebenen 
Trajcktorientypen sind gut zu erkennen. Die Trennlinien zwischen den Bereichen entsprechen den 
Gleichuneen 15.261 (I— II), 15.311 (I— III) und 15.371 (II— III) und stimmen sehr gut mit dem gefunde- 
nen qualitativen Verhalten iiberein. Der Ubergang zwischen dem Platzen durch Kalte (II) und 
Uberdehnung (III) wird nur durch die Modellierung als scharfe Linie erkennbar. In der Praxis 
fafit man diesen Bereich besser als einen stetigen Ubergang auf. Zur Orientierung sind die Kar- 
tenballons aus der Versuchsreihe „Glaschke" (siehe Tab. 16. 1| ) mit eingetragen. Sie zeigen deutlich, 
dafi fiir die verwendeten Ballons die maximale Flugweite von der Kalteemphndlichkeit des Ballon- 
gummis abhangt, wahrend die PlatzgroBe zweitrangig ist. Dieser Befund laBt sich auch auf die bei 
Ballonwettbewerben verwendeten Modell^l iibertragen. 

1 Als Standardwert wird in dieser Arbeit g = 9, 81 m/s 2 verwendet. 
2 Ballongr6Ben iiber 25 cm sind bereits ausreichcnd. 
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Die Flugweite bei fester Leermasse und verschiedenen Startauftrieben gibt einen tieferen Ein- 
blick in die Abhangigkeit der Flugweite von der Ballonfullung (Abb. [521 • Dcutlich crkcnnbar sind 
die bereits angesprochenen drei ausgezeichneten Flugweiten d max , dburst und <iburst2- Die jeweils 
grofite erreichte Flugweite wird durch d max markiert. Links davon befinden sich Trajektorien mit 
Ballons, die nach dem Erreichen der Scheitelhohe intakt absteigen. Unmittelbar nach rechts fallt 
die Flugweite auf ciburst ab - hier fallen die Ballons geplatzt zu Boden. Bei noch grofierem Auf- 
trieb nimmt die Scheitelhohe und Flugweite weiter ab (c?burst2)- Besonders deutlich wird die starke 
Abhangigkeit der maximalen Flugweite von der Ballonfullung. Eine Genauigkeit von mindestens 
0,5 1 beim Fiillen ist notwendig, und selbst dann kann die Wetterabhangigkeit von ^opt die ge- 
naue Lage um bis zu einem Gramm verschieben. Der Verlauf der Scheitelhohen kann aus einem 
Schnitt mit Mq = 7, 1 g in Abbildung 15.11 entnommen werden. 

Die numerische Losung der Differentialgleichungen bietet auch die Moglichkeit, den EinfluB dcs 
vernachlassigten Drucksprungs auf die Ballontrajektorie zu diskutieren. Zum Vergleich wurde ein 
optimal gefiillter Heliumballon mit und ohne Drucksprung gerechnet sowie ein leicht iiberfulltcr 
Ballon fAbb. [5~5]) . Der endliche Drucksprung Ap macht sich in einer ausgepragten Asymmetrie des 
Ballonaufstiegs bemerkbar: Beim Aufsteigen wird der Ballon starker komprimiert (relativ zu einem 
Ballon mit Ap = 0), so dafi der Auftrieb schneller abnimmt und der Ballon regelrecht gebremst 
wird. Beim anschliefienden Abstieg wird dieser „gespeicherte" Auftrieb wieder freigesetzt und 
verlangsamt den Sinkfiug. Dieser Effckt wird zusatzlich durch die starke Temperaturabhangigkeit 
der Permeabilitat verstarkt[f|. Dieser Mechanismus erklart auch, warum die Flugweite ciburst weit 
geringer ist als die Halfte von d max . Trotz dieser starken Vereinfachung der analytischen Rechnung 
behalten die abgeleiteten Skalierungsgesetze ihre Giiltigkeit (vgl. Abb. 15.4( 1. 

AbschlieBend wird noch der optimale Auftrieb fur verschiedene Leermassen in Abbildung 15.41 
dargestellt. Leermassen kleiner 10 g entsprechen Fall I und sind fur das Erreichen grofier Flugweiten 
besonders zu empfehlen. Anhand der numerischen Kurven kann auch die Proportionalitatskonstan- 
te in Gl. 15.261 bestimmt werden: 




(5.40) 



Als Faustformel gilt, dafi der Auftrieb rund 20 % der Leermasse betragen sollte. 



3 Roberts 1821 geht von einer anderen Temperaturabhangigkeit der Permeabilitat aus und findet nur eine schwach 
ausgepragte Asymmetrie. 
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Trajektorien Heliumballons 
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Abbildung 5.3: Trajektorien fur einen optimal gefullten Ballon, einen leicht iibervollen Ballon und 
einen optimal gefullten Ballon ohne Drucksprung. 
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Kapitel 6 

Ballonfluge 



6.1 Einleitung 

Der Entwurf ernes Trajektorienmodells mufi mit Experimenten verglichen werden, um die verwen- 
deten Annahmen beurteilen zu konnen und Hinweise auf unbekannte Phanomene zu erhalten. Die 
praktische Durchfiihrung von Experimenten wird durch mehrere Faktoren erschwert: Kartcnbal- 
lons haben nur einen geringen Auftrieb, der es unmoglich macht Messtechnik anzubringen, ohne 
die Ballonparameter zu beeinfiussen. Dies schliefit detaillierte Informationen iiber die Flugbahn 
aus. Werden vor dem Start eines Ballons alle Parameter erfafit, sind die einzigen zuganglichen 
Informationen die Balloneigenschaften beim Start und der Landeplatz. Dieser Ansatz wird jedoch 
durch die geringe Fundquote von im besten Fall 10 % (vgl. Tab. 16.11) erschwert, die mehrere tausend 
Ballons erfordert, um statistisch signifikante Aussagen treffen zu konnen. 

Aus diesem Grund wurde zunachst eine Vorstudie unter kontrollierten Bedingungen durch- 
gefuhrt (siehe Kapitel I6.3[) und durch Daten von Ballonwettbewerben erganzt. Der Vorteil von 
Ballonwettbewerben besteht in den umfangreichen Daten, die oft iiber mehrere Jahre mit meh- 
reren Dutzenden gefundenen Karten pro Wettbewerb anzugeben sind. Der Nachteil besteht in 
den unbekannten Ballonparametern, zumal die wissenschaftliche Auswertung erst im Nachhinein 
crfolgt und typischerweise nicht Teil der Wettbewerbsplanung ist. Eine besondere Datcnqucllc 
besteht in der Zusammenarbeit mit Herrn Haberlandt [37], der Ballonstarts iiber Jahre privat 
durchgefuhrt hat und genaue Angaben zu den Startbedingungen machen kann. 

Die Eckdaten jedes gefundenen Kartenballons bestehen somit aus Startort und Startzeit sowie 
dem Fundort. Einblicke in den Ballonflug und die beteiligte Physik konnen nur indirekt gewon- 
nen werden, indem jeder Fundort mit einer individuell berechneten Trajektorie verglichen wird. 
Unstimmigkeiten geben dann einen Hinweis auf vernachlassigte Effekte oder Fehler des Trajekto- 
rienmodells. 

Die Anzahl von mehreren tausend Funden macht eine automatische Auswertung zwingend er- 
forderlich. Zunachst wurden alle verfiigbaren Daten einheitlich formatiert und die entsprechenden 
Trajektorien automatisch berechnet. Die Berechnung und Auswertung der Trajektorien erfolgte 
mit einem selbstcntwickcltcn Computcrprogram 

6.2 Funddaten 

Alle verwendctcn Daten sind in Tabelle IQ1 zusammengefafit. Angegeben ist jeweils der Zeitraum, 
in dem die Ballonstarts erfolgt sind, die Anzahl der Starttage in dieser Zeit sowie die Gesamtzahl 
der gefundenen Ballons. In den letzten Spalten ist das verwendete Fiillgas und die Fundquote 
angegeben. 

1 Anhang [B] stellt das entwickelte Programm kurz vor. 
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Referenz 


Zeitraum 


Tage 


Ballons 


Gas 


Fundquote 


Quelle 


Campingkirche 


2000 - 


2008 


46 


881 


He 


10-15 % 


m 


Glaschke 


1998 - 


2001 


7 


7 


H 2 


11 ± 4 % 




Haberlandt A 


2003 - 


10.2004 


60 


162 


He 


4,2 ±0,3 % 


m\ 


Haberlandt B 


10.2004 - 


08.2007 


215 


640 


He 


3,6 ±0,1 % 


m\ 


Haberlandt C 


09.2007 - 


2008 


72 


281 


He 


4,3 ±0,3 % 




Trostberg 


1998 - 


2010 


13 


240 


He 


9,2 ±0,6 % 


[53 



Tabelle 6.1: Auflistung aller verwendeten Ballonwettbewerbe mit Angaben zu den gefundenen 
Ballonkarten. Geschatzte Fundquoten sind ohne statistischen Fehler angegeben. 

„ Campingkirche" und „Trostberg" sind regelmafiig stattfindende Ballonnugwettbewerbe, die 
auch aktuell noch durchgefiihrt werderd. Sie bilden iiber Jahre einen homogenen Datensatz. Die 
Daten „ Trostberg" enthalten auch Informationen zu der Zeit, die zwischen dem Start eines Ballons 
und dem Fund der Karte vergangen ist. 

Die umfangreichste Datenquelle ist „ Haberlandt". Die Vielfalt der durchgefuhrten Starts mach- 
te eine Aufteilung nach den Angaben des Experimentators in drei BlScke notig: Block A enthalt 
explorative Versuche mit BallonvoluminaH bis zu 17 1, die schrittweise auf etwa 8,2 1 reduziert 
wurden. Als Mittelwert werden 13 1 angenommen. Die Ballons wurden iiber den Tag verteilt bis 
in den friihen Abend gestartet. Block B sind Starts mit der als optimal erkannten Fiillung von 

8.2 1 und in Block C wurden Wettervorhersagen des Windes und des Niederschlags mit einbezogen, 
um die Flugweite weiter zu erhohen und die Verlustrate durch ungunstigcs Wetter zu minimieren. 
Starts erfolgten hier nur noch am friihen Abend. 

„ Glaschke" reprasentiert eigene Ballonexperimente, in denen neben der Uhrzeit, dem Gewicht 
und Auftrieb auch die Ballons selbst photographisch erfafit wurden. 

6.3 Eigene Experimente 



Nr 


Datum 


Startzeit 


Falg} 


Mofe] 


"z,o[km/h] 


d[km] 


^Flug 


1 


19.09.1998 


18:25 


7,1 


5,4 


5,0 


181 


? 


2 


03.08.1999 


14:23 


8,6 


7,5 


5,4 


30 


< 2 h 


3 


16.08.1999 


14:28 


5,8 


7,2 


4,6 


105 


? 


4 


18.09.1999 


19:45 


5,8 


7,4 


4,8 


143 


? 


5 


21.09.1999 


16:13 


5,5 


7,2 


4,6 


136 


? 


6 


01.01.2000 


14:46 


6,4 


7,2 


4,8 


21 


34 min 


7 


28.01.2001 


12:35 


4,3 


7,2 


4,2 


267 


? 



Tabelle 6.2: Startdaten der zuriickgeschickten Ballonkarten. 

Tabelle IST21 enthalt die Startdaten aller gefundenen Kartenballons. Die letzten beiden Spalten 
enthalten die Flugweite und Informationen iiber die Flugzeit, soweit sie aus den Angaben des 
Finders ableitbar waren. Die Ballons erreichten teilweise eine Flugweite von mchr als 250 km, 
wobei die Weite mit abnehmendem Auftrieb des Ballons zunimmt. 

Die bekannten Parameter der Ballons erlauben es, eine vollstandige Simulation jeder Ballon- 
trajektorie durchzufiihren. Als einziger freier Parameter wurde die Hohe, in der der Ballon durch 
unbekannte Einniisse platzt, verwendet. Mit Flugweiten um 100 km bewegen sich die Trajektori- 
en noch unter der Auflosungsgrenze des Windmodells (siehe Kap. [3]), so dafi die freie Platzhohe 

2 Daten nach 2008 waren teilweise erst nach dem Beginn der Auswertung verfugbar und sind dicser Arbeit noch 
nicht berticksichtigt. 

3 Es sind jeweils die geometrischen Volumina angegeben. 
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Abbildung 6.1: Vergleich der Flugweiten der Ballons aus Tabelle RT21 mit Simulationen des Refe- 
renzballons. 

solange verandert wurde, bis die Modellflugweite dem Experiment entsprichl|f|. Die abgeleiteten 
Platzhohen sind zusammen mit der Temperatur in dieser Hohe und der Gesamtflugzeit in Tabelle 
16.31 enthalten. Die Platzhohen fallen grob in zwei Gruppen: Ballons, die in rund 10 km Hohe bei 
Temperaturen um — 50°C mit einem Volumen von etwa 30 1 platzen und Ballons, die deutlich 
friiher platzen. Die erste Gruppc cntspricht den Erwartungen an die Haltbarkeit des Ballongum- 
mis (Kap. [2]), so dafi keine weiteren aufieren Einniisse angenommen werden miissen. Ballons der 
zweiten Gruppe platzen deutlich friiher, typischerweise bei wenigen Minusgraden, was den Einflufi 
von Eis auf der Ballonhiille naheleg1@. 

Die Simulationen zeigen aber auch weitergehende Unterschiede in den beiden Fallen mit be- 
kannter Flugzeit: 

• Bei Ballon 6 konnte der Finder die Landung der Karte beobachten, so dafi die Flugzeit genau 
bekannt ist. Aus der Flugzeit von 34 Minuten folgt eine maximale Hohe von etwa 1,2 Kilo- 
metern, was deutlich unter der berechneten Hohe ist. Die langere Flugzeit des berechneten 
Ballons deutet auf eine Unterschatzung der Winde durch die interpolierten Radiosondenda- 
ten hin. 

• Die Flugzeit von Ballon 2 weicht leicht von der gemessenen Flugzeit ab, die in diesem Fall 
aber nur eine Obergrenze darstellt (wann der Ballon genau vor dem Fund gelandet ist, ist 
nicht bekannt). Im Rahmen der erwarteten Genauigkeit ist der Fehler akzeptabel. 

Abbildung 16.11 gibt eine abschliefiende Ubersicht aller Flugweiten zusammen mit einer Simulati- 
onsreihe des Standardballons aus Tabelle 15.11 Hier wird besonders deutlich, dafi alle gestarteten 
Ballons weit vom berechneten Optimum entfernt sind. Der Wert dieser Versuchsreihe besteht vor 
allem in der Absicherung der Bedingungen, unter denen einen Ballon platzt. 



4 Den Fundort exakt anzupassen ist im Rahmen der Modellgenauigkeit nicht moglich. 
5 Dieses Phanomen ist bei Radiosondenballons bekannt I77| . 
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Nr 




7min[°C] 


*Flug|>] 




1 


10,0 


-47 


4,0 


27,7 


2 


5,8 


-14 


2,3 


24,0 


3 


5,3 


-14 


2,3 


19,0 


4 


10,8 


-54 


4,4 


31,0 


5 


7,6 


-32 


3,2 


22,1 


6 


4,0 


-12 


1,7 


17,1 


7 


9,8 


-53 


4,4 


27,5 



Tabelle 6.3: Auswertung der zuriickgeschickten Ballonkarten. 



6.4 Startdynamik 

Fiir den erfolgreichen Start eines Ballons ist das Windprofil in Bodennahe von groBem Ein- 
flufi. Wahrend der Wind Ballons urn geschlossene Hindernisse fast storungsfrei mitnimmt, stellen 
durchlassige Objekte wie Baume und Antennen mogliche Hindernisse nach dem Start dar. 

Zur genaueren Beurteilung wurde ein Teil der in dem vorhergehenden Kapitel geschilderten Bal- 
lonstarts verwendet, urn das bodennahe Windprofil zu bestimmen. Die Vorgehensweise entspricht 
der Windmessung durch Pilotballone: Nach dem Start wird der Ballon in regelmafiigen Abstanden 
angepeilt und der H5hen- und Azimutwinkel notiert. Aus der bckannten Steigrate v z des Ballons 
kann dann die dreidimensionale Flugbahn und somit die Windgeschwindigkeit bestimmt werden. 
Als Modellfunktion wird das bereits vorgestellte logarithmische Windprofil (Gl. I3.5[) verwendet, 
da keine Quantifizicrung der atmospharischen Stabilitat durchgefuhrt werden konntc: 

v(z) = v Q ln(z/z ) (6.1) 

Zur kompaktcrcn Darstcllung wurden alle Vorfaktorcn zu cincr Skalcngcschwindigkcit vq zusam- 
mcngezogen. Mit Hilfe dieses Windprofils kann die Flugbahn des Ballons im Raum als Funktion 
der Hohe z explizit bestimmt werden: 

r(z) = ^(z\n(z/z a )-(z-z )) (6.2) 

r ist der horizontal gemessene Abstand zum Startort. Fiir die Abschatzung der Kollisionsgefahr 
mit umstehenden Objekten ist der Abflugwinkel <fi besser geeignet: 

tan(^) = -Lr (6.3) 
r(z) 



1 



hx(z/z ) -lv 



(6.4) 



ps 20°— (6.5) 

Zur Bestimmung der Skalengeschwindigkeit vq wurden fiir mehrere Ballonstarts Geschwindig- 
keitsprofile erstellt ( Abb. und mit Gl. l6.1l verglichen. Da eine einzelne Trajektorie keine gemit- 
telte, reprasentative Messung darstellt, ist ein logarithmisches Windprofil nicht fiir alle Messungen 
klar erkennbar. Dies trifft besonders auf Starts mit boigem Wind zu. Aus diesem Grund konnte 
nur ein Parameterbereich ermittelt werden: 

v(z) - (l,5...5,0)xln(-^) ^ (6.6) 

\2 m/ h 

Der Wertebereich deckt sich mit unabhangigcn Bcstimmungcn des Windprofils durch Radioson- 
den |69] . Der zu erwartende kleine Abflugwinkel von optimal gefiillten Kartenballons macht eine 
besondere Sorgfalt beim Start erforderlich. 
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6.5 Ballonwettbewerbe 

Die Auswertung der weiteren Daten aus Tabelle UTTl erfolgt auf eine andere Weise. Da die genau- 
en Parameter der gestarteten Ballons (insbesondere das Ballongewicht, der Startauftrieb und in 
vielen Fallen die Uhrzeit) nicht bekannt sind, ist eine detaillierte Auswertung von jedem Start 
nicht moglich. Statt dessen wird der folgende Ansatz gewahlt: Fiir alle Starts wird pauschal ein 
Kartenballon mit der Masse 5,1 g und einer Karte von 2 g angenommen. Unbekannte Start- 
zeiten werden durch den Schatzwert 15:00 Uhr ersetzt. Der unbekannte Startauftrieb wird durch 
einen Ensembleansatz umgangen: Fiir jeden Ballonstart wird eine ganze Reihe virtueller Fluge mit 
F a ,o = 1 ■ ■ ■ 4 g durchgcfiihrt und die jeweils am besten passende Trajektorie zur Fehlerberechnung 
verwendet. Dieses Vorgehen fiihrt im Prinzip zu einer Uberschatzung der Trajektoriengenauig- 
keit, aber in der Mehrzahl der Falle ist der Einflufi des Auftriebs auf die Flugrichtung so gering, 
daB der wahre Fehler kaum beeinflufit wird. Zusatzlich wird auch fiir jeden Start der optimale 
Auftrieb berechnet, der die Flugweite maximiert. Die folgcnden zwei Abschnitte gehen naher auf 
die Flugweiten und Modellfehler ein. Abschliefiend werden starke Abweichungen der Trajektorien 
diskutiert. 



6.6 Ballonflugweite 

Die individuellen Flugweiten d der zur Verfiigung stehenden Kartenballons bilden die Grundlage, 
um die in den vorangegangenen Kapiteln beschriebenen Effekten zu untersuchen und quantitativ zu 
erfassen. Hierbei muB das Augenmerk auf einer sorgfaltigen Zusammenfiihrung der verschiedenen 
Flugtage und Wettbewerbe liegen, da die pro Flugtag verfugbaren Daten nur in den seltensten 
Fallen ausreichen, um eine aussagekraftige Statistik zu erstellen. 

Die Schwierigkeit besteht in den individuell verschiedenen Wetterbedingungen jeden Flugtags, 
die einen einfachen Vergleich der Daten ausschliessen. Am prominentesten ist die zeitliche Verande- 
rung des Windfelds, welche direkt die Flugweite eines Kartenballons beeinflufit. Die Eigenschaften 
des Windfelds erlauben aber eine elegante Losung dieses Problems. Zum einen ist die Mehrzahl 
der Trajektorien nur schwach gekriimmt - dies zeigen berechnete Trajektorien und Kartenbal- 
lonfundorte - so daB die Betrachtung auf das vertikale Windprofil v(z) beschrankt werden kann. 
Zum anderen wird eben dieses Windprofil im Mittel gut durch eine lineare Funktion beschrieben 
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Verteilung Flugweiten, Wildon, 1999-2008 




Abbildung 6.3: Verteilung der charakteristischen Flugweiten d max und (iburst- 




(vgl. Abb. 13. lj) . Dies lafit erwarten, dafj sich der Einflufi des Windes durch eine einfache Skalierung 
der Flugweite korrigieren lafit. Urn diesen Ansatz weiter zu stiitzen, wurden die Trajektorienrech- 
nungen aus Kapitel 16.101 verwendet, um die Verteilung der charakteristischen Flugweiten c? max 
und (iburst zweidimensional darzustellen. Trotz der unterschiedlichen Trajektorien, die mit den 
jeweiligen charakteristischen Flugweiten verkniipft sind (vgl. Abb. I5.3[) . zeigt die Verteilung in 
Abbildung 16.31 eine sehr starke Korrelation. 

Dicse Argumente erlauben die Einfuhrung cincr skalierten Flugweite d, die aus den individuell 
berechneten Flugweiten c? max und (iburst fur jeden Kartenballon berechnet wird. Die stiickweise 
linear skalierte Entfernung d ist definiert durch: 

falls d < dburst 

d := { falls d>d max (6.7) 

sonst. 

Die Wahl des kritischen Verhaltnisses (iburst /d max = 0,3 wurde uber ein Mittelwert aus vielen 
Simulationen (vgl. Abb. 16.31) bestimmt. Die Skalierung der Flugweite dburst auf einen festen Wert 
wird durch die Aussagekraft dieses Wertes motiviert: Flugweiten groBer als dburst fordern eindeutig, 
dafi der Ballon den Scheitelpunkt der Trajektorie intakt durchflogen hat, wahrend kleinere Werte 
wahrscheinlich (aber nicht notwendigerweise) auf einen geplatzten Ballon hindeuten. Normicrte 
Flugweiten grofier Eins sollten im Sinne dieser Definition nicht vorkommen. Eine genaue Analyse 
dieses Falls erfolgt mit einer ausfiihrlichen Fehlerbetrachtung in den Kapiteln 16.81 und 16.91 

Die Verwendung einer skalierten Flugweite erlaubt, windbedingte Unterschiede in den Flugwei- 
ten zu korrigieren, was die Kombination von Kartenballondaten verschiedener Flugtage ermoglicht. 
Dariiber hinaus ist die Mittelung verschiedener Flugtage eine Moglichkeit, Aussagen unabhangig 
und ohne genaue Kenntnis des Wettergeschehens der einzelnen Starttage treffen zu konnen. 

Die im Folgenden vorgestellten Verteilungen der skalierten Flugweite erlauben auch die Um- 
gehung eines weiteren Problems: Wenn der Startauftrieb eines Ballons nicht bekannt ist (z.B. im 
Falle der Ballonwettbewerbe) , kann aus der Flugweite eines einzelnen Ballons nicht eindeutig auf 
die Flugzeit geschlossen werden (vgl. Abb. 15.21) . Bei der Betrachtung einer Verteilung kann aber 
jeder skalierten Flugweite d im Mittcl cine zuverlassigere Flugzeit iFiug zugeordnet werden: 



tpiue ~ 20 h x 



0,45 x (J/0,3) 2 / 3 falls d < 0,3 

-0,E 
0,7 



0,45 + 0,55^^ sonst 
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Verteilung Flugweiten 
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Abbildung 6.4: Kumulative Verteilungen der normierten Flugweite. 

Diese Gleichung wurde ebenfalls aus einem geeigneten Mittel iiber viele Simulationen gewonnen. 

Zunachst wird die Verteilung der normierten Flugweiten in Abbildung 16.41 untersucht. Der 
Idealfall, d. h. alle Ballons erreichen die maximal mogliche Flugweite, ist gegeben durch 

P (>dj = e(i-J) (6.9) 

mit der Heaviside-Sprungfunktion O. Eine erste Inspektion der Verteilungen zeigt bereits, dafi 
dieses idealisierte Szenario nicht einmal annahernd eine treffende Beschreibung gibt. Zur weite- 
ren Untersuchung wurden die Verteilungen anhand ihrer Ahnlichkeit in zwei Gruppen aufgeteilt: 
„H a berlandt B" und „C" (Gruppe 2) sowie die verbleibenden Datensatze (Gruppe 1). Die mei- 
sten Flugweiten in Gruppe 1 liegen deutlich unter ciburst, so dafi hier in Ubereinstimmung mit 
den Angaben zu „Haberlandt A" von sehr stark gefiillten Ballons auszugehen ist. Uberraschcnd 
ist die starke Streuung der Flugweiten, die sich nicht alleine mit unterschiedlich gefiillten Ballons 
crklarcn lafit. Gerade im Bereich stark gefullter Ballons variiert die Flugweite nur schwach mit 
dem Auftrieb. 

Gruppe 2 zeigt deutlich grofiere normierte Flugweiten, was durch die gezielte Optimierung des 
Startauftriebs zu erwarten war. Aber welche Verteilung ist in diesem Fall iiberhaupt zu erwarten? 
Das Volumen der Ballons in Gruppe 2 bewegt sich im Rahmen des technisch moglichen dicht am 
optimalen Wert 8,2 1. Dennoch ist es unwahrscheinlich, dafi alle Ballons die maximale Flugweite 
erreichen, da der optimale Bereich sehr schmal ist und die genaue Lage vom Wettergeschehen 
abhangig ist. In diesem Licht ist es sinnvoll, eine annahernde Gleichverteilung der Ballonflugweiten 
zu postulieren: 

p(>dj = l-d (6.10) 

Aber selbst diese vielleicht etwas zu pessimistische Verteilung iiberschatzt den Anteil grofier ska- 
lierter Flugweiten. 

Es miissen also weitere Einflufifaktoren vorliegen, die bisher vernachlafiigt wurden und weit- 
gehend unabhangig von den Balloneigenschaften operieren. Der universelle, annahernd exponen- 
tielle Verlauf aller Verteilungen deutet in die gleiche Richtung. Bereits genannte Kandidaten sind 

39 







Haberlandt A - 






Haberlandt B 


* 1 




Haberlandt C 




***** *'*"•, 


Campingkirche 






Trostberg 




d burst 


i i i i 




UV/Ozon (Kap.[H3J und Regen fKap. 1572"]) . Zur genaueren Analyse dieser Einflufifaktoren wurde 
fiir die Daten „Haberlandt A" (hier sind die Startzeiten bekannt) jeweils eine eigene Statistik fur 
die Starts vor 16:00 Uhr und danach erstellt. Um einen quantitativen Vergleich der Verteilungen 
in Gruppe 1 mit einem idealisierten Szenario zu ermoglichen, wurde zusatzlich eine Jahresstatistik 
fiir den Startort Wildorj^ berechnet. An jedem Tag des Jahres 2004 wurde ein 13 1 Heliumbal- 
lon (vgl. Tab. 15. ip virtuell um 14:00 Uhr(UT) gestartet. Abbildung 16.51 stellt noch einmal alle 
Verteilungen aus Gruppe 1 und die berechnete Jahresstatistik dar. Obwohl ein konstantes Bal- 
lonvolumen die in der Praxis auftrctcndcn Variationen nicht beriicksichtigt, ist die Diskrepanz zu 
den empirisch bestimmten Verteilungen signifikant. Alleine „Haberlandt A>16:00" zeigt einen mit 
der Jahresstatistik vergleichbaren Einbruch in der Verteilung bei d = 0, 25 durch zu stark gefiillte 
Ballons. 

Es ist hilfreich, den bereits geschilderten Befund fiir eine genauere Diskussion zu quantifizieren. 
Dazu wird der Ballonflug in zwei Prozesse zerlegt: Den ungestorten (d. h. idealisierten) Ballonflug, 
der alleine durch die Balloneigenschaften und die Windstromung bestimmt wird, und einen von 
aufien eingreifenden Mechanismus, der den Ballon zerstort bzw. zu Boden bringt. Unter der Annah- 
me, dafi beide Prozesse unabhangig voneinander sind, k5nnen die jeweiligen Uberlebensfunktionen 
cinfach miteinander multipliziert werden: 

F(d) = F T rajoct(d) X F Decay (d) (6.11) 

Im hier betrachteten Fall entspricht die Gesamtiiberlebensfunktion F der Verteilung der skalierten 
Entfernungen p(> d), und die Teilfunktion i*T> a ject bezieht sich auf den ungestorten Ballonflug. 
-Fbecay beschreibt die Zerstorung des Ballons durch die bereits angesprochenen Umwcltcinfliisse. 
Die treffende Modellierung dieser Effekte erscheint auf den ersten Blick unmoglich, da keine de- 
tailliertcn Informationcn iiber die zugrundclicgcndcn Mcchanismcn vorliegen. Durch die Mittclung 
verschiedener Flugtage und die chaotische Natur des Wettergeschehens ist es aber sinnvoll anzu- 
nehmen, dafi diese Prozesse unabhangig und homogen in Zeit und Ort operieren. In diesem Fall 

"Startort der Daten „Haberlandt" 
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Referenz 


r[h] 


Kommentar 


Trostberg 


1,5-2,5 


Starts Nachmittags, Juni, Juli 


Habcrlandt A < 16:00 


3- 4 


Starts ab Morgens 


Campingkirche 


4- 5 


Starts Nachmittags 


Glaschke 


4- 6 


Starts Nachmittags 


Habcrlandt A > 16:00 


10- 12 


Starts bis in den Abend 


Habcrlandt B 


12- 16 


Starts Abends 


Haberlandt C 


20- 26 


Starts Abends, Wetterbericht 



Tabcllc 6.4: Ermittcltc Lcbcnsdauern fur verschiedene Datcnsatzc. 



crhalt man zwangslaufig einen exponentiellen Zcrfall 

-Fbccay(d) = exp(-t Flug (J)/T) (6.12) 



wobei die Flugzeit ^Fiug gemafi Gl. 16.81 eingesetzt wird. Alles relevante Wissen iiber die externen 
Effekte ist in der Lebenszeit r enthalten. Es ist klar, dafi Gl. 16.121 nur eine erste Naherung ist, 
die durchaus vorhandene Korrelationen wie z. B. den Tagesrhythmus aufier Acht lafit. Durch die 
bereits vorhandenen Unsicherheiten der Daten ist aber eine weitergehende Interpretation nicht 
sinnvoll. 

Das Einsetzen von Gl. 16.121 in Gl. I6.11l zur Bcstimmung einzelner Lebenszeiten r erfordert die 
angemessene Wahl einer ungestorten Uberlebensfunktion -Prrajcct- Hierfiir mufi notwendigerweise 
auf Simulationen zuruckgegriffen werden, die durch ihre Konstruktion keine externen Einfliisse 
beinhalten. Fiir Gruppe 1 wird die bereits genannte Jahresstatistik 2004 verwendet, wahrend fur 
die deutlich verschiedene Gruppe 2 Gl. 16.101 eingesetzt wird. Die Durchfuhrung einer regelrechten 
Ausgleichsrechnung zur Bestimmung der Lebenszeiten wiirde den bereits genannten Limitierung 
nicht gerecht werden, so dafi statt dessen jeweils von Hand mogliche Wertebereiche ermittelt 
wurden. 

In Tabelle l5?4l sind alle ermittelten Lebenszeiten zusammengestellt. Im Zweifel wurden grofiziigi- 
gere Wertebereiche angegeben, um vor allem die Unsicherheit von -Frraject besser abzubilden. Durch 
die genau bekannten Ballonparameter des Datensatzes „ Glaschke" konnte die Lebenszeit direkt 
aus Tabelle 16.31 ermittelt werden. Trotz der unterschiedlichen Datcnqucllen und Methoden zeigen 
„Haberlandt A<16:00", „ Campingkirche" und „Glaschke" mit vergleichbaren Startbedingungen 
eine bemerkenswerte Ubereinstimmung der Lebenszeit von 4-5 Stunden. „Haberlandt A> 16:00" 
sowie „B" mit Starts am spaten Nachmittag bzw. friihen Abend sind ebenfalls in zufriedenstellen- 
der Ubereinstimmung, vor allem wenn man den deutlichen Unterschied in den Ballonfiillungen und 
Flugzeiten bedenkt. Mit dem Datensatz „Haberlandt C", der zusatzlich nur Starts bei trockenem 
Wetter enthalt, wird die grofite Lebenszeit erreicht. Mit „Trostberg" ist die kiirzeste Lebenszeit 
im Stundenbereich gegeben. 

Abschliefiend werden die Angaben in Tab elle [6 . 4 1 verwendet . um jeweils den einzelnen Effekten 
Lebensdauern zuzuordnen. In erster Naherung ergeben sich die aufgefuhrten Lebensdauern aus 
einer Summe der verschiedenen Beitrage: 

1 « -^ + ^ + ^- (6.13) 

T TLicht TRcgcn TRest 

Aus dem Vergleich der beiden „Haberlandt A"-Bl6cke ergibt sich der durch Tageslicht bestimmte 
Anteil zu TLicht ~ 10 h. Vergleiche mit „Haberlandt B" und „C" erlauben die Zuordnung fiir 
Regen T Kcgen « 20 h und die verbleibenden (d. h. unbekannten) Effekte TR CSt ~ 20 h. Solange diese 
Aufteilung nicht durch eine direkte Verkniipfung der Ballonflugweiten mit Wetterdaten bestatigt 
wird (vgl. Kap. l3T2|) . sind diese Angaben vorlaufige Schatzwerte. 

In Anbetracht dieser Aufteilung bleibt alleine der Datensatz „ Trostberg" ratselhaft, der eine 
deutlich kleinere Lebensdauer aufzeigt. Nach Rucksprache mit dem Veranstalter wurde die Aus- 
wertung noch cinmal mit cincm vcranderten Ballonmodell wiederholt. Das Ergebnis anderte sich 
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Referenz 


% Intakt 


Campingkirche 


7,8 


Glaschke 


0,0 


Haberlandt A 


4,8 


Haberlandt B 


37,4 


Haberlandt C 


56,3 


Trostberg 


0,0 



Tabelle 6.5: Anteil intakt gelandeter Ballons. 

aber nur geringfiigig, ohne diesen Befund erklaren zu konnen. Moglicherweise fiihrten die Start- 
zeiten in den Monaten Juni und Juli zu einem besonders starken Einflufi der Sonneneinstrahlung, 
aber auch die Nahe des Startorts zu den Alpen ist ein Umstand, der von dem Trajektorienmodell 
nicht detailliert beriicksichtigt wird. 

6.7 Intakte Ballons 

Die kumulativen Verteilungen in Abbildung 16.41 erlauben auch eine Abschatzung des Anteils der 
Ballons, die intakt gclandet sind. Ballons, die weiter als (iburst geflogen sind, miissen den Scheitel- 
punkt intakt passiert haben. Da die Hiille beim Abstieg wieder schrumpft und dadurch zunehmend 
entlastet wird, ist es sehr wahrscheinlich, dafi ein Grofiteil dieser Ballons ebenfalls intakt landet. 
Der aus den kumulativen Verteilungen entnommene Anteil dieser Ballons ist in Tabelle 16.51 als 
Schatzwert zusammengestellt. 

6.8 Fehler 

Das vorgestellte Modell muB durch eine Fehlerbetrachtung gepriift und beurteilt werden. Un- 
stimmigkeiten zwischen dem Fundort eines Kartenballons und der korrespondierenden Trajektorie 
konnen zwei prinzipiell verschiedene Ursachen haben: Zum einen kann das verwendete Trajektori- 
enmodell selbst fehlerbehaftet sein, und zum anderen konnen die mit dem gefundenen Kartenballon 
verkniipften Daten ungenau sein. 

Der Fundort eines Kartenballons wird selten prazise angegeben. In der Regel wird die nachste 
Stadt genannt, was im Mittel eine Genauigkeit von wenigen Kilometern erlaubt (vgl. Kap. [4]). 
Wenn ein Kartenballon jedoch ins Meer oder in ein ffiefiendes Gewasser fallt, hat auch ein ge- 
nauer Fundort nur noch eine begrenzte Aussagekraft: Durch den zusatzliche Transport geht jede 
Information iiber den wahren Landeplatz verloren. In diesem Fall sind fast beliebig grofie Abwei- 
chungen moglich sind. Oft ist auch die genaue Startuhrzeit eines Kartenballons nicht bekannt. Hicr 
mufi mit Schatzwerten gearbeitet werden, die aber um bis zu mehrere Stunden falsch sein konnen. 
Besonders bei dynamischen Wettersituationen ist dann ein deutlicher Einflufi auf die berechnete 
Trajektorie zu erwarten. Zuletzt besteht die Moglichkeit, dafi beim Aufzeichnen der Daten Fehler 
aller Art auftreten. Das konnen Verwechslungen und Zahlendreher bei der Archivierung sein, Pro- 
blems bei der Identifizierung des Fundorts oder gezielte Falschungen durch den Finder oder den 
Teilnehmer eines Ballonwettbewerbs. 

Das Trajektorienmodell besteht aus mehreren Komponenten, die unterschiedlich stark die Ge- 
nauigkeit ciner Trajektorie bccinflussen. Die wichtigsten Teile sind die zugrunde liegenden Glci- 
chungen (Kap. [5]), die verwendeten Wetterdaten (Kap. [3]) und die numerische Integration des 
Gesamtsy stems. Von zentraler Bedeutung sind Fehler in den Grundgleichungen. Diese konnen je- 
doch nur indirekt erschlossen werden, indem gepriift wird, ob der Gesamtfehler alleine durch die 
begrenzte Genauigkeit der Wetterdaten erklart werden kann. Das numerische Losungsverfahren 
und der Zeitschritt sind bereits so gewahlt, dafi der Integrationsfehler vernachlassigbar ist. 
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Abbildung 6.6: Das durch die Trajektorie mit minimalem A m i n aufgespannte Koordinatensystem. 
Der offene Kreis ist der Endpunkt der Trajektorie. Der korrespondierende Fundort bcfindct sich 
im Abstand d zum Startort. 

Die begrenzte Genauigkeit des Wettermodells ist nicht alleine auf fehlerhafte Einzelmessungen 
zuriickzufuhren, sondern vor allem auf die geringe Abtastung der Wetterwerte in Zeit und Ort |92] . 
Windbewegungen variieren bereits im Minuten- und Stundenbereich |97] , was deutlich kiirzer als 
das tibliche Startintervall der Radiosonden von 12 Stunden (bzw. in seltenen Fallen 6 Stunden) ist. 
Diese Unterabtastung verursacht zwei Fehlerbeitrage: Zum einen konnen schnell variierende Wind- 
beitrage gar nicht durch das interpolierte Wettermodell wiedergegeben werden (Auflosungs-Effekt), 
zum anderen werden schnelle Fluktuationen durch die Unterabtastung in niedrigere Frequenzbe- 
reiche verschoben (Alias-Effekt) und andern so ebenfalls das Modell. Wahrend die hochfrequenten 
Bcitrage im Laufe einer Trajektorienberechnung gemittelt und damit gedampft werden, summieren 
sich langperiodischen Bcitrage systematisch auf. 

Die gleiche Argumentation kann auf die Abtastung im Ort angewcndct werden. Bei einem 
mittleren Abstand der verwendeten Radiosondcnstationen von bis zu 200 km und typischen Kor- 
relationslangen des Windfelds von mehreren 100 bis 1000 km [57] ist der EinfluB aber etwas geringer 
ausgepragt. Zu beachten ist, daB sich bei kurzen Trajektorien die kleinskaligen Windstrukturen 
nicht herausmitteln konnen und damit den Fehler dominieren. 

Dartiber hinaus wird die erreichbare Genauigkeit durch fehlende Wetterdaten fundamental 
limitiert: Interpolationsverfahren konnen die Auswirkungen fehlender Daten mildern, aber keines- 
falls aufheben. In diesem Licht sind eventuell durch das Interpolationsverfahren selbst eingefuhrte 
Fehler zu vernachlassigen. 

Der Gesamtfehler A einer Trajektorie als Funktion der Flugzeit t ergibt sich aus der Summe 
des Auflosungsfehlers Ar os und des Alias- Fehlers AAiias^ 

Alt) = A Rcs (t) + A AVms (t) (6.14) 

Bcdingt durch den systematischen Charakter wachst AAiias ebenso wie die Flugweite d naherungs- 
weise linear mit der Zeit: 

A A lia S « t(Av) (6.15) 

d{t) w t(v) (6.16) 

(Av) ist der mittlere Geschwindigkeitsfchler und (v) die mittlere Windgeschwindigkeit. Da die 
Ursache des Windfehlers (Av) eine Verschiebung im Zeit- bzw. Ortsspektrum ist, ist (Aw) oc (v) 
zu erwarten. Fiir die meisten Trajektorien liefert Aft es nur einen kleinen Beitrag, weshalb auf die 
genaue Analyse des zeitliche Verlaufs verzichtet wird. 

Zur Fehlerauswertung wird der Fundort mit der Vorhersage der Simulation iiber ein speziel- 
les Koordinatensystem in Verbindung gebracht (siehe Abb. 16. 6p . Die y-Koordinate wird aus der 
Lange des Lots auf die Trajektorie bestimmt (positiv fiir eine Abweichung nach links vom Star- 
tort aus gesehen), wobei die Trajektorie an den Enden iiber ihre Tangente verlangert wird. Die 
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Abbildung 6.7: Gleitendes Mittel des absoluten Trajektorienfehlers A m ; n . Gl. 16.201 ist als Gerade 
mit eingetragen. 



Lange der Trajektorie bis zum Fufipunkt ist die in Flugrichtung positiv gezahlte x-Koordinate. 
Ein einfaches und robustes Fehlermafi ist die kiirzeste Entfernung A m j n des Kartenfundorts x / zur 
berechneten Trajektorie, die in der Abbildung durch eine gestrichelte Linien dargestcllt ist. Aus 
den vorangegangenen Ausfuhrungcn erscheint der relative Fehler des Fundorts besonders geeignet, 
der ebenfalls berechnet wird: 



mm{\x f 

Aniii-, 



£ Trajektorie} 



^Rcl ■ — 



(6.17) 
(6.18) 



Durch die Definition von A m ; n ist der relative Fehler beschrankt: Ar c i < 1. 

Nun werden die in Gl. 16.171 und Gl. 16.181 angegebenen Fehlermafie berechnet und statistisch 
ausgewertet. Zuerst wird die Abhangigkeit von A m i n von der Flugweite d untersucht. Um stati- 
stische Schwankungen zu reduzieren, werden die einzelnen Fehlerwerte A m j ni fe aufsteigend sortiert 
und ein gleitendes Mittel angewendet: 



1 



2/ 



k=-l 



A 



min,i+fc 



(6.19) 



Der Mittelungsbereich umfafit 41 Wcrtc (/ = 20) und wird an den Randern symmetrisch verkiirzt. 
Abbildung 16.71 fafit alle Fehler zusammen. Trotz der unterschiedlichen Datenquellen zeigt sich ein 
universelles lineares Ansteigen des absoluten Fehlers mit der Flugweite. Dies bestatigt die be- 
reits ausgefiihrte Uberlegung, dafi der Alias-Effckt in der Windfeldabtastung einen systematischen 
Fehler bedingt. Das Fehlerverhalten lafit sich gut durch 



A min (d) = 0,13 d 



(6.20) 



beschreiben. Der mittlcrc Fehler von 13 % pafit gut zu den Ergebnissen von [541, IM] . wobei zu 
beriicksichtigen ist, dafi sich Gl. 16.201 nur auf einen eindimensionalen Fehler senkrecht zur Trajek- 
torie bezieht. 
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Abbildung 6.8: Relativer Fehler der einzelnen Trajektorienrechnungen. Gl. l6.21l ist mit eingetragen. 



Da der Alias-Effekt von dem Anteil hoher Frequenzen in den Windfluktuationen abhangt, ist es 
durchaus moglich, dafi Wettersituationen vorliegen, die auch mit einer Abtastung von 12 h relativ 
gut zur rekonstruieren sind. Um einen Einblick in die Fehlerverteilung zu erhalten, wird der relative 
Fehler nach Gl. 16.181 fur jeden Kartenballon ausgewertet und in einer kumulierten Verteilung 
aufgetragen. Trajektorien kiirzer als 200 km wurden ausgeschlossen, um einer Verzerrung der 
Verteilungen durch den vernachlassigten Auflosungsfehler vorzubeugen. 

Wie die Schwankungen im mittleren Fehler erwarten lassen, sind hier die Ergebnisse weit 
uneinheitlicher (Abb. Qualitativ folgen aber alle Verteilungen annahernd der exponentiellen 
Verteilung 

F(A Rcl ) = exp(-^ (6 .2i) 

die mit dem Mittelwert in Gl. 16.201 konsistent ist. Die Abnahme der verfiigbaren Daten mit zu- 
nehmender Flugweite erschwert die Beurteilung, ob diese Verteilung wirklich unabhangig von der 
Flugweite d ist. Zumindest im Bereich d < 1000 km unterstiitzen die Daten diese Annahme. 

Da kerne zeitlich hoher aufgelosten Winddaten zur Verfugung stehen, laBt sich dieses Fehlerver- 
halten nicht eindeutig mit dem Windspektrum in Verbindung bringen. Einzelne Messungen und 
theoretische Arbeiten deuten aber in diese Richtung [57J HI] • 



6.9 Ausreifier 

Ballons, die deutlich weiter fliegen als berechnet oder deren Trajektorie einen grofien Fehler auf- 
zeigt, sind mogliche Hinweise auf Schwachen oder Fehler in dem vorgestellten Trajektorienmodell. 
Um zu entscheiden, ob es sich in diesen Fallen nur um Extremwerte der beschriebenen Fehlerver- 
teilung (Abb. I6.8[) oder wirkliche Probleme handelt, wurden alle auffalligen Ballons noch einmal 
genauer untersucht. Ballons mit einer normierten Flugweite d iiber 1,1 wurden in einer Grup- 
pe zusammengefafit. Aus den verbleibenden Ballons wurden alle Ballons mit einem Fehler A m i„ 
iiber 200 km in einer weiteren Gruppe zusammengefafit. Um der Ursache der Abweichungen auf 
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Referenz 


Ballons 


d> 1,1 


1) 


2) 


3) 


4) 


A min > 200 km 


1) 


2) 


3) 


4) 


Campingkirche 


881 


2 


1 


1 








6 





1 


4 


1 


Glaschke 


7 










Haberlandt A 


156 






















Haberlandt B 


636 


13 


4 


9 








31 


7 


17 


3 


4 


Haberlandt C 


266 


7 


7 











19 





15 





4 


Trostberg 


211 











Tabelle 6.6: Auffallige Daten. Die detaillierte Beschreibung der Kategorien 1-4 befindet sich auf 
Seitelli 

die Spur zu kommen, wurden fur jeden Ballon noch einmal Trajektorien „von Hand" berechnet 
und die Struktur des interpolierten Windfelds untersucht. Anhand dieser Ergebnisse wurden vier 
Fehlerkategorien erstellt: 

1. Vertretbare Abweichung. Bei Gruppen sind die Ergebnisse fur die restlichen Ballons gut oder 
die Flugweite war nahe 200 km. 

2. Schwierige Windfeldrekonstruktion (Divergenzen, Scherstromungen, keine dominante Haupt- 
stromung). Dies umfafit sowohl Orts- als auch Zeitintcrpolationsprobleme. 

3. Starke Abweichung von der berechneten Trajektorie und anderen Ballons (falls verftigbar). 
Flug gegen die vorherrschende Windrichtung ohne erkennbaren Grund. 

4. Verbleibende starke Abweichungen. 

Tabelle 16.61 enthalt die Klassifizierung aller auffalligen Ballons. Der Grofiteil der Abweichungen 
ist dabei der ungenauen Rekonstruktion des Windfelds anzulasten (Fall [5]). Fiir die Gruppe der 
Ballons mit normierten Flugweiten uber 1,1 deckt diese Erklarung alle Fliige ab, so daB aus den 
Daten kein Widerspruch zur berechneten maximalen Flugweite abgeleitet werden kann. 

Es verbleibt aber ein Teil Ballonfiuge mit unklaren Ursachen (Fall [4]). In einigen Fallen konnten 
die Ballons an Kusten angespiilt worden sein, aber auch durch starke Regenfalle o. a. veranderte 
Trajektorien konnten in diesen Fallen eine Rolle spielen. 

Besonders ratselhaft sind die Falle, in denen Ballons gegen die vorherrschende Windrichtung 
gedriftet sind, ohne daB die Windfeldrekonstruktion problcmatisch crschcint (Fall[3]). Hicr konnten 
auch einfache DatenfehleiQ, fchlcrhaft ausgefiillte Ballonkarten oder Verwechslungen bei der Be- 
stimmung des Fundorts vorliegen. Bei Ballonwettbewerben bleibt schliefilich auch die Moglichkeit 
der gezielten Falschung. Herr Haberlandt nennt eine Falschung (oder einen Scherz) durch den 
Finder [sic] der Karte. Aber selbst wenn man Fall [3] durchgangig als gezielte Falschungen betrach- 
tet, bleibt der auf die Teilnehmer der Ballonwettbewerbe hochgerechnete Anteil sehr gering. Auch 
wenn eine abschliefiende Beurteilung zu diesem Zeitpunkt nicht moglich ist, stellen diese Falle mit 
einem Anteil um 0,5 % nur einen Bruchteil aller Kartenballons dar. 

6.10 Jahresstatistiken 

Die absolute Flugweite eines Kartenballons wurde bis jetzt kaum diskutiert, da sie in ho hem 
MaBe von den Wind- und Wetterbedingungen abhangt und nur wenig uber die Giiltigkeit eines 
Trajektorienmodells aussagt. Diese fiir die Ausrichter von Ballonwettbewerben und naturlich deren 
Teilnehmer wichtige GroBe soli nun naher beleuchtet werden. 

Mit dem erfolgreich validierten Modell wurden „virtuelle" Ballontrajektorien und Landepunkte 
berechnet, die eine genaue Vorstellung von den erreichbaren Flugweiten ermoglichen. Dieser Ansatz 

7 Bei der Formatierung der Ballonflugdaten wurden je nach Datcnqucllc in 0,1 % bis 1 % der Falle Fehler erkannt 
(Zahlendreher, namensgleiche Fundortc u. a.). 
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Verteilung Flugweiten, Wildon, 1999-2008 
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Abbildung 6.9: Berechnete Jahresstatistik. Zum Vergleich zeigt S(d) die Verteilung in Gl. 16.251 



erlaubt eine weit bessere Statistik, die mit echten Ballons kaum mit vertretbarem Aufwand moglich 
ware. 

Als Startpunkt dieser virtuellen Fliige dient Wildorjf], das auch der Startort des grofiten Da- 
tensatzes „Haberlandt A-C" ist. Starts wurden fiir jeden Tag der Jahre 1999 - 2008 mit der 
Startzeit 14:00 Uhr(UT) durchgefiihrt. Bei einer Flugzeit von rund 20 Stunden liegt der Zeitab- 
lauf der Trajektorie symmetrisch zu den Wetterdaten um 00:00 Uhr(UT), wodurch systematische 
Einniisse der Zeitinterpolation reduziert werden. Da das Atmospharenmodell deterministisch ist, 
crgibt cin Ballon pro Startzeitpunkt jeweils eine eindeutige Trajektorie. Fiir jeden Startzeitpunkt 
wurde der Ballonauftrieb iterativ verandert, um die charakteristischen Flugweiten c£ max und dbm-st 
automatisch zu bestimmen. 

Die zweidimensionale Verteilung der Fundorte mit maximaler Flugweite d max auf der Erdober- 
flache wird gut durch eine bivariate Gaufi- Verteilung beschrieben. Bedingt durch die Westwinddrift 
liegt der Modalwert 660 km weiter ostlich bei (46°15'N 24°4'E). Trotz dieser Asymmetrie lassen 
sich die in Abbildung 16.91 zusammengefassten kumulierten Abstandsverteilungen gut durch einfa- 
che Rayleigh-Verteilungen annahern: 

<r mH = 911 km (6.22) 

Cburst = 258 km (6.23) 

Die Form der Abstandsverteilung geht im wesentlichen auf die zugrundeliegenden rayleighverteiltcn 
Windgescliwindigkeiten zuriick. Bei den grofiten Abstanden weichen die Daten etwas von dieser 
einfachen Verteilung ab, aber dieser Bereich wird bereits durch starke statistischen Schwankungen 
dominiert. Die gewonnenen Verteilungen spiegeln bis jetzt nur den Idealfall wieder - jeder Ballon 
erreicht nach dem Start die maximal mogliche Flugweite. Mit einer bekannten Uberlebensfunktion 
_F(siehe Kap. l6.6l) . bezogen auf die normierte Flugweite, kann diese Verteilung aber naherungsweise 

s Wildon/Steicrmark, 46°54'N 15°33'E 



S'max(> d) = exp 

<Sburst(> d) = exp 
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Verteilung Flugweiten 
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Abbildung 6.10: Erreichte Flugweiten. Gl. 16.251 dient hier ebenfalls als Vergleichspunkt. 



korrigiert werden, urn ein realistischeres Bild zu erhalten: 

S(>d) = f S{d/x)F{x)dx (6.24) 



Exemplarisch sei hier das Ergebnis fiir den Fall gleichverteilter normierter Flugweiten gegeben: 
S max (>d) = cxp (-^-) - \/|^- erfc f -J^— ) ( 6 ' 25 ) 

V 2cr max/ V2cr max V V 2<7 max / 



o 2 



max 



-=-5- +0 (6.26) 



Bcrcits dieses cbenso idcalisierte Szenario zeigt eine deutlich geringere Wahrscheinlichkeit grofier 
Flugweiten. 

Der Vergleich der Kartenballondaten mit dieser Betrachtung zeigt, dafi alle Ballonwettbewerbe 
weit hinter der Modellfunktion Gl. 16.251 zuriickbleiben. Nur Haberlandt B und besonders Haber- 
landt C erreichen beachtliche Flugweiten. Hier ist aber zu beachten, dafi gezielt Tage mit hohen 
Windgeschwindigkeiten zum Start ausgewahlt wurden, wahrend in der berechneten Jahresstatistik 
iiber alle Tage des Jahres gleich gemittelt wurde. 

Diese Statistik schliefit immer wieder gemeldete extreme Flugweiten wie Fliige iiber den Atlan- 
tik mit Sicherheit aus. Alte Ausgaben des Guinness-Buch der Rekorde verzeichnen noch „Rekord- 
nugweiten" deutlich iiber 10 000 km [55], aber in neueren Aunagen ist diese Kategorie folgerichtig 
nicht mehr enthalten5. Auf der Basis der hier durchgefiihrten Rechnungen liegt die maximal er- 
rcichbarc Distanz (zumindest in Europa) wenig iiber 3000 km. 



9 Diese Kategorie wird mindestens seit 1995 nicht mehr gefuhrt, Mitteilung der Redaktion des Guinness World 
Records Buches vom 23.2.2010. 
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6.11 Zeit bis zum Fund 



Abschliefiend wird noch die Zeit betrachtet die vergeht, bevor eine Karte gefunden und zuriick- 
geschickt wird. Unter der Annahme, dafi die Fundwahrscheinlichkeit an einem Landeplatz zeit- 
unabhangig ist, erhalt man eine exponentielle Verteilung des Anteils der noch nicht gefundenen 
Karten Ff unc i: 

F lnnd (t) = exp(-t/r) (6.27) 

Die mittlere Liegezeit r bis zum Fund der Karte betragt fur die Daten von Glaschke 6,7 Tage, 
nach den Daten Trostberg ergibt sich eine Lebensdauer von 6-18 Tagen. Werden nur kurze Zeiten 
im Wochenbereich betrachtet, liefert Gl. 16.271 eine sehr gute Beschreibung, auch wenn signifikan- 
te Abweichungen erkennbar sind. Besonders deutlich wird das bei den Daten von Haberlandt, 
die Karten beinhalten, die nach mehr als einem Jahr zuruckgeschickt wurden. Hier erfordert die 
Heterogenitat der Fundorte, von einer ganzen Verteilung von mittleren Liegezeiten auszugehen: 

F lund {t) = / cxp(-t/r)/(r)dr (6.28) 
Jo 

Im konkreten Fall der Daten von Haberlandt sind bereits zwei Liegezeiten ausreichend, um eine 
gute Naherung zu erhalten: 

F fund (*) w o,85fflq>r-^J + 0,15 6X13 ( _ 22^d) (6 ' 29) 

Es ist zu erwarten, dafi die Fundwahrscheinlichkeit mit der Bevolkerungsdichte und der Siedlungs- 
struktur in Zusammcnhang gebracht werden kann. Diese Idee wurde aber aufgrund der unzurei- 
chenden Datenlage nicht weiter verfolgt. 
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Kapitel 7 

Diskussion 



Das in dieser Arbeit entwickelte Modell ermoglicht vielfaltige Einblicke in die Dynamik eines Kar- 
tenballons. Der Vergleich der durchgeffihrten Simulationen mit Ergebnissen von Ballonwettbewer- 
ben erlaubt fundierte Aussagen fiber die typische Ballontrajektorie: Ballons erreichen regelmafiig 
Hohen iiber 12 km, widerstehen dabei Tempcraturen bis zu —60 °C und konnen bis zu 24 Stun- 
den und langer in der Luft bleiben. Als wichtigster Parameter bei der Flugweitenmaximierung 
hat sich die Ballonfiillung herausgestellt. Der berechnete optimale Startauftrieb von rund 20 % 
der Leermasse ist dabei in guter Ubereinstimmung mit dem unabhangig experimentell ermittelten 
optimalen Ballonvolumen von 8,2 1 [37]. Ist die Ballonfiillung optimiert, haben alle anderen Ballo- 
neigenschaften bemerkenswerterweise nur noch einen geringen EinfluB auf die Flugweite. Alleine 
die maximal erreichbare Flughohe zeichnet sich durch einen fast quadratischen Einflufi aus. Re- 
gen und UV-Strahlung (besonders im Zeitraum Juni - Juli) erscheinen als die starksten externen 
Einfliisse, die die Lebenszeit eines Ballons auf wenige Stunden beschranken konnen. 

In der Praxis sind unter diesen Umstanden Flugweiten bis zu 3000 km moglich, was die Angaben 
in [55] deutlich nach oben korrigiert. Nach der Kenntnis des Autors liegt der in dieser Arbeit 
durch Trajektorienrechnungen bestatigte Rekord bei 2321 km [37]. Der auf den ersten Blick gering 
crscheinende Drucksprung eines Ballon tragt besonders zu einer langen Verweilzeit in groficn Hohen 
und damit zu grofien Flugweiten bei. Leider bleiben Ballonwettbewerbc durch zu stark gefiilltc 
Ballons und Starts tagsiiber im Sommer weit unter diesen Werten. 

Erfreulich sind die im Vergleich zu den Betrachtungen in Kapitel U sehr hohen Fundquoten 
(siehe Tab. 16. ll) . die vermuten lassen, dafi die Mehrheit der gefundenen Karten auch anschliefiend 
zuriickgeschickt wird. 

Trotz der erfolgreichen Validierung des Kartenballonmodells besteht noch Raum fur Verbes- 
serungen. Die Qualitatskontrollc und Integration der primaren Radiosondendaten konnte starker 
systematisiert werden [35]. Prinzipiell ist auch die Einbindung eines digitalen Gelandemodells in 
das Atmospharenmodell denkbar (z. B. [3B]). Ein solcher Schritt erfordert aber zwangslaufig ei- 
ne konsistente Integration des Windfelds, was weit iiber die in dieser Arbeit verwendete einfache 
Interpolation hinausgehen wiirde. 

Alle Aussagen iiber das Ballonvcrhaltcn im Flug sind nur indirekt erschlossen. Dirckte Ein- 
blicke wiirde nur eine Flugsimulation in einer Klimakammer bringen (in eingeschrankter Form 
von |52j verwendet), die auch definitive Aussagen fiber das Materialverhalten von Gummi bei tie- 
fen Temperaturen und starker Dchnung liefern wfirden. Ebenso waren Laborexperimente notig, 
um das Zusammenspiel von UV-Strahlung und Oxidantien wie Ozon bei der Wirkung auf den 
Ballongummi zu untersuchen. 

Zur Steigerung der Zuverlassigkeit ist auch die Verwendung von Ballongespannen moglich, 
die weitaus vielfaltigere Moglichkeiten zur Optimierung bieten. Aber auch bei einer weitgehenden 
Optimierung eines Kartenballons bleibt ein Problem bestehen: Ein normaler Kartenballon verffigt 
fiber keine Regelmechanismen, die das Aufsteigen fiber die maximal mogliche Hohe und damit das 
Platzen aktiv verhindern. 
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Anhang A 

Nutzliche Konstanten 



9 


9,81 m/s 2 


Erdbeschleunigung 


-RErde 


6371 km 


Mittlerer Erdradius 


e 


1,602176 487(40) x 10~ 19 C 


Elemcntarladung 


k B 


1,380 6504(24) x 10~ 23 J/K 


Boltzmann Konstante 


R 


8,314 472(15) J mor 1 FT 1 


Gaskonstante 


a 


5,670 400(40) x lO" 8 ^ 


Stefan-Boltzmann-Konstante 


c 


299 792 458 m/s 


Lichtgeschwindigkcit 


h 


6,626 068 96(33) x 10~ 34 Js 


Plancksches Wirkungsquantum 


N A 


6, 022 141 79(30) x 10 23 mor 1 


Avogadro Konstante 


v m 


24,464 0424 1 mor 1 


Molares Volumen bei Standardbedingungen 


T 


-273,15°C 


Absoluter Nullpunkt 


NTP 


25°C 101325 Pa 


Standardbedingungen (nicht genormt) 


STP 


0°C 101325 Pa 


Normbedingungen 



Materialeigenschaften 





p[kg/m 3 ] 


%[Luft] 




AS 




Luft 


1,204 


100,000 


28,964 


0,0240 


1,005 


He 


0,1785 




4,0026 


0,15 


5,19412 


H 2 


0,0837 




2,016 


0,1861 


14,3 


N 2 


1,166 


78,084 


28,013 




1,040 


2 


1,332 


20,942 


31,999 




0,9196 


co 2 


1,842 


0,038 


44,010 




0,8504 


Gummi 


1000? 






0,16? 


1,4? 



Dichten jeweils bei 20°C NP. Kinematische Viskositfit von Luft : 1,5x10 5 m 2 /s 
Spezifische Dichte verschiedener Materialien: 



PE-Folic 


0,12 


g/dm 2 


Holzstab 2 mm 


0,44 


g/dm 


Al-Folie (Schokoladenverpackung) 


0,30 


g/dm 2 


Al-Folie 


0,37 


g/dm 2 


Strohhalm 


0,11 


g/dm 


Papier 


0,83 


g/dm 2 


Buttcrbrotpapicr 


0,43 


g/dm 2 



52 



Anhang B 



Datenverarbeitung 



am gearbeiten Extras 



Laden | 28.11.2008 |Schwape Karle _l]|lB _d Zumck | Vqi | 3h _^J|5QQQ 




|fertig 36° 53' 6" [J 16° 49' 50" O 



Benutzeroberflache des selbst entwickelten Auswertungsprogramms Balloon Analyser. Alle inner- 
halb des Programms referenzierten Objekte sind links in einem Baum iibersichtlich zusammenge- 
stcllt. Aufgefiihrt werden externe Wetterdaten und Kartenmaterial (Ordner Umwelt), importicrtc 
Dateien mit Ballonfunden (Ordner Ballonfunde), sclbstdcfinicrtc Startorte fur unabhangige Si- 
mulationen (Ordner Orte), benutzerdcfinicrtc Ballons zur Durchfiihrung virtueller Fliige (Ordner 
Ballons), Simulationen zum Erstellen von Statistiken und Fundauswertungen (Ordner Simula- 
tion) und schliefilich Analysen zum automatischen Erstellen verschiedener Diagramme (Ordner 
Analysen). 

Die Objekte im Baum konnen zu einem Ablaufschema verkniipft werden, das automatisch 
alle Simulationen berechnet und auswertet. Neuberechnungen mit veranderten Parametern oder 
aktualisierten Daten sind dadurch bequem moglich. 

Das Abspeichern aller Einstellungen in einer Projektdatei dokumentiert die aufgesetzten Si- 
mulationen und ermoglicht das spatere Nachvollziehen der durchgefiihrten Berechnungen. 
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Anhang C 

Formelzeichen 



V : Volumen 

A : Oberflache 

A q : Qucrschnittsflache 

R : Radius / Gaskonstantc 

Rq : Radius druckfreicr Ballon 

D : 2R Durchmesser/Diffusivitat 

e : Volumcnkorrckturfaktor 

VJnax : Maximalcs Ballonvolumen 

VGummi : Gummivolumcn dcr Ballonhiille 

Mcummi: Gewicht der Ballonhiille 

p : Dichte 

N : Stoffmenge [mol] 

T : Temperatur 

t : Zcit 

T : Strahlungsflufi 

A : Warmclcitfahigkcit / Geographische Lange 

/3 : Geographische Breite 

Nu : Nusselt-Zahl 

Re : vL/v Reynolds-Zahl 

Pr : ijCp/X Prandtl-Zahl 

Gr : Grashof-Zahl 

v : Kinematischc Viskositat 

j] : Dynamische Viskositat 

c w : Widerstandsbeiwert 

F a : Nettoauftrieb (= Auftrieb — Gewichtskraft) 

Q : Permeabilitat (bezogen auf Partialdruck) 

6 : Permeabilitat (bezogen auf Stoffmenge) 

kff : Henry-Konstante 

n : Teilchendichte 

p(> x) : Komplementare Verteilungsfunktion 

m : Molekularmasse 

A : Fehler 

v : Gcschwindigkeit 

r : Abstand 
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Anhang D 

Reibungskoeffizient 




0.01 1— hi— i 1— kL_i 

10" 2 10" 1 10° 10 1 10 2 10 3 10 4 10 5 10 6 



10' 



Re 



Formel von Weidenschilling 



24/Re fur Rc < 1 

24 Re" ' 6 fur 1< Rc < 800 Re = 

0,44 fur 800 < Re 



2Rv 



(D.l) 



Fit dcr Datcn von 123 : 



24 

— + 7, 6 Re' 
Rc 



0,693 



0.44 



(D.2) 
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Anhang E 

Mittleres Atmospharenprofil 



Hohc[km] 


T[°C] 


u[m/s] 


Druck[hPa] 


Dichte [p/po] 


0,02 


8,6 


5 


1013,8 


1,00 


0,50 


6,8 


8 


956,5 


0,95 


1,00 


4,7 


9 


899,1 


0,90 


2,00 


0,0 


10 


795,1 


0,81 


3,00 


-5,1 


12 


700,4 


0,73 


4,00 


-10,7 


14 


616,3 


0,65 


5,00 


-16,8 


17 


540,3 


0,59 


6,00 


-23,2 


20 


473,6 


0,53 


7,00 


-30,0 


23 


410,9 


0,47 


8,00 


-37,0 


26 


354,8 


0,42 


9,00 


-44,2 


29 


306,8 


0,37 


10,00 


-51,6 


32 


264,1 


0,33 



Tabelle E.l: Mittlere Daten der Atmosphare als Funktion der Hohe nach [T3] . 
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